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INTRODUCTION 


Use  of  nondestructive  evaluation  techniques  for  assessing 
criticality  of  defects  in  aircraft  structural  components  is  an 
attractive  method  for  defining  their  residual  strength  or  life¬ 
time.  It  also  provides  a  basis  for  defining  regular  inspection 
intervals  for  a  structural  element  during  its  service  and  for 
determining  the  need  for  mandatory  repair.  Application  of  such 
techniques  to  components  made  of  fiber  reinforced  composites  is 
being  studied  by  various  investigators.  A  number  of  nondes¬ 
tructive  inspection  techniques  have  been  and  are  being  examined 
for  detecting  flaws  in  composite  structural  members  and  quanti¬ 
fying  their  location,  size,  and  shape.  These  include  ultra¬ 
sonics,  acoustic  emission,  x-ray  and  thermography,  use  of  pene¬ 
trants  as  well  as  holography. 

After  the  flaw  characteristics  have  been  quantified  by 
means  of  ND  measurements,  it  is  necessary  to  assess  criticality 
by  the  use  of  analytical  techniques.  For  metal  structural  com¬ 
ponents  such  techniques  are  well  known  and  their  validity  is 
well  established.  Since  fiber  composites  have  been  in  use  for 
only  a  few  years,  suitable  analytical  methodologies  for  the 
definition  of  critical  flaws  are  not  fully  developed.  Also,  be¬ 
cause  of  their  highly  anisotropic  properties,  stress  analysis  in 
a  flawed  composite  component  involves  lengthy  and  complicated 
computations  and  different  types  of  flaws  may  result  in  a  wide 
variety  of  failure  modes.  For  these  reasons  usefulness  of 
analytical  methods  for  assessing  flaw  criticality  in  composite 
structures  has  not  been  clearly  demonstrated. 

Since  it  is  practically  impossible  to  use  a  single  analyti¬ 
cal  technique  to  cover  the  whole  spectrum  of  fatigue,  fracture 
and/or  failure  for  various  types  of  flaws,  a  logical  procedure 
for  studying  criticality  is  to  select  a  few  types  of  commonly 
encountered  flaws  and  develop  appropriate  analytical  models  for 
each  type  of  defect.  Such  analytical  studies  should,  of  course, 


be  supplemented  with  nondestructive  as  well  as  destructive  tests 
to  examine  the  range  of  validity  of  the  model.  Such  attempts 
have  been  made  in  references  1  and  2  for  interlaminar  disbonds  in 
a  laminated  beam.  These  defects  are  quite  common  in  laminated 
members  because  they  can  exist  as  "birth"  defects  or  can  be 
created  during  service  by  foreign  object  impact  or  various  other 
reasons.  Presence  of  a  disbond  does  not  affect  the  structural 
performance  of  a  laminated  beam  or  plate  member  under  in-plane 
loads  except  when  such  loads  can  cause  local  buckling  due  to  in¬ 
plane  compression  or  shear.  Under  transverse  shear, catastrophic 
and  slow  growth  of  a  disbond  may  occur  due  to  quasi-static  and 
cyclic  loading,  respectively.  Criticality  and  growth  rate  of 
disbonds  in  laminated  beams  under  quasi-static  or  cyclic  trans¬ 
verse  shear  have  been  determined  (ref.  2)  by  the  use  of  princi¬ 
ples  and  methods  of  fracture  mechanics  as  well  as  empirical 
crack  growth  laws.  Failure  of  a  disbonded  laminate  under  com¬ 
pressive  loading  has  been  defined  in  terms  of  buckling  and  an 
elastic  stability  analysis  has  been  utilized  for  assessing  criti¬ 
cality.  Sizes  and  shapes  of  disbonds  implanted  in  test  specimens 
have  been  quantified  and  their  growth  monitored  by  the  use  of 
ultrasonic  "C"  scans.  Results  from  tests  have  been  compared  with 
analytical  predictions  and  usefulness  of  the  NDE  methodology  has 
been  demonstrated.  An  approach,  which  is  qualitatively  similar 
to  the  one  described  above,  has  been  utilized  in  reference  3  for 
studying  growth  of  inherent  disbonds  ("birth"  defects)  in  lamin¬ 
ated  composites  subjected  to  compression  fatigue.  Although  the 
philosophy  and  the  analytical  methodology  of  assessing  criticality 
of  disbonds  have  been  shown  to  be  effective,  their  range  of  valid¬ 
ity  as  well  as  their  usefulness  for  more  complicated  disbond  geom¬ 
etries  and  loading  conditions  are  yet  to  be  demonstrated.  The 
work  reported  in  this  study  is  an  attempt  in  this  direction. 

In  the  present  work,  ultrasonic  "C"  scans  are  used  to  quanti¬ 
fy  shapes  of  disbonds  with  nearly  elliptic  planforms  which  are  im¬ 
planted  in  wide  beams.  "C"  scans  are  also  used  to  monitor  growth 


of  the  disbonds  under  cyclic  transverse  shear.  Destructive  tests 
are  performed  to  determine  quasi-static  failure.  Sandwich  beams 
with  compression  skins  containing  circular  as  well  as  nearly  el¬ 
liptic  disbonds  are  loaded  to  failure  for  determining  criticality 
of  the  disbonded  skin. 

Analytical  methods  are  developed  for  determining  the  effect 
of  elliptic  shaped  disbonds  on  the  stress  field  in  a  composite 
laminate.  The  analysis  yields  the  strain  energy  release  rates, 

Gj-r  Git»  and  Gjjj  »  associated  with  opening  (Mode  I),  sliding 
(Mode  II),  and  tearing  (Mode  III)  modes,  respectively,  as  func¬ 
tions  of  position  along  the  boundary  of  the  disbond.  The  possi¬ 
bility  of  predicting  flaw  growth  under  cyclic  loading  is  examined. 
Effects  of  in-plane  compressive  or  shear  loads  are  also  included 
in  the  analysis,  and  therefore,  failure  of  disbonded  laminates 
due  to  elastic  instability  is  also  predicted.  Data  correlation 
studies  are  performed  with  the  help  of  experimental  data. 

In  addition,  wave  propagation  studies  are  performed  on  resin 
and  laminate  samples,  for  quantifying  changes  in  wave  velocity 
and  attenuation  induced  by  moisture  conditioning  and  fatigue 
damage.  Results  are  studied  to  examine  the  usefulness  of  velo¬ 
city  and  attenuation  measurements  as  NDT  techniques  for  quantita¬ 
tive  assessment  of  such  damages. 


ANALYTICAL  METHODS 


DESCRIPTION  OF  THE  MODEL 

A  laminated  composite  plate  of  finite  thickness  with  a  de¬ 
lamination  between  any  two  of  the  laminae,  as  shown  in  figure  1, 
is  considered  in  this  study.  The  displacements  and  stresses  in 
such  a  plate,  due  to  applied  loads,  can  be  calculated  as  the  sum 
of  the  solutions  of  two  problems  as  described  below. 

1.  The  problem  of  calculation  of  stresses  and  displace¬ 
ments  without  the  presence  of  the  flaw  -  self-equili¬ 
brating  surface  tractions  on  the  two  surfaces  of  the 
delamination,  as  determined  from  this  solution,  are 
used  as  inputs  for  the  second  problem. 

2.  Effects  of  tractions  equal  in  magnitude  but  opposite  in 
sign  to  those  calculated  in  (1) ,  acting  on  the  surfaces 
of  the  delamination. 

In  many  cases  the  solution  to  the  first  problem  is  either 
known  or  can  be  determined  by  known  analytical  techniques  using 
two-dimensional  elasticity  or  laminated  plate  theories.  When 
the  periphery  of  the  delamination  is  located  quite  far  from  ti.e 
edges  of  the  plate,  the  solution  to  the  second  problem  is  insen¬ 
sitive  to  the  boundary  conditions  at  these  edges  and  therefore, 
the  flaw  is  assumed  to  be  located  in  a  plate  of  finite  thickness, 
but  the  other  two  dimensions  of  the  plate  are  infinitely  large. 

A  rigorous  three-dimensional  elasticity  solution  of  the 
second  problem  under  consideration  seems  possible,  but  will  be 
extremely  complicated  because  of  the  anisotropy  of  the  individual 
layers  as  well  as  the  laminated  structure  of  the  plate.  The  com¬ 
plexity  of  the  problem  will  be  clear  from  past  studies  on  a  two- 
dimensional  elasticity  solution  for  the  problem  of  a  laminated 
composite  beam  containing  delaminations  (refs.  2,4)  and  a  three- 
dimensional  elasticity  solution  for  an  elliptic  crack  in  an 


anisotropic  medium  (refs.  5,6).  In  this  study  a  lower  order 
structural  theory,  based  on  a  laminated  plate  theory  with  shear 
deformations  (ref.  7)  which  is  modified  as  described  in  Appendix 
A- 1  to  include  the  effects  of  tractions  on  the  bounding  surfaces 
of  the  plate  (see  fig.  2),  is  used  to  obtain  results  for  the 
following  cases. 

1.  Contributions  of  mode  I,  mode  II  and  mode  III  types  of 
deformations  to  strain  energy  release  rates  as  func¬ 
tions  of  position  on  the  periphery  of  a  disbond  with 
elliptic  plan  form  as  shown  in  figure  1. 

2.  Critical  values  of  membrane  stresses  (compressive,  as 
shown  in  fig.  1,  or  shear)  to  cause  buckling  failure 
of  the  unsupported  part  of  the  compression  skin  above 
the  disbond. 

Under  the  assumptions  of  shear  deformation  plate  theory  and 
Griffith  separation,  the  analysis  for  case  (1)  is  self-consistent 
and  yields  an  exact  solution.  The  model  consists  of  two  lamin¬ 
ated  plates,  as  shown  in  figure  1,  bonded  to  each  other  every¬ 
where  except  over  the  area  of  the  disbond.  As  described  in 
Appendix  A-2,  the  theory  admits  the  presence  of  three  distinct 
singularities  at  the  periphery  of  the  disbond  associated  with 
normal  and  two  tangential  stresses  in  the  form  of  line  loads  due 
to  discontinuities  in  displacement  gradients  across  the  disbond 
boundary,  thus  yielding  finite  values  of  G^,  G.^  and  Gjjj,  the 
strain  energy  release  rates  corresponding  to  modes  I,  II  and  III, 
respectively.  These  types  of  singularities  have  been  noted  in 
the  dynamic  problem  of  debonding  of  a  Timoshenko-beam  from  a 
rigid  substrate  (ref.  8).  Quasi-static  growth  of  disbonds  have 
also  been  studied  by  several  authors  (refs.  8-18)  using  various 
types  of  structural  theories.  Results  from  references  18  and  1 
are  compared  in  references  2  and  «  with  a  2-D  elasticity  solution 
for  a  laminated  beam  containing  delaminations  in  an  environment 
dominated  by  transverse  shear  stresses.  The  comparison  shows 
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that  the  strain  energy  release  rates  obtained  from  the  two  theor¬ 
ies  are  close  to  each  other  when  the  ratio  of  disbond  length  to 
the  thickness  of  each  of  the  two  beams  in  the  disbonded  laminate 
is  greater  than  five.  For  decreasing  values  of  disbond  length, 
the  lower  order  theory  results  become  non-conservative.  This 
phenomenon  is  expected  because  of  inherent  limitations  in  the 
shear  deformation  beam  theory.  It  is  noted,  however,  that  for 
small  values  of  disbond  length,  flexural  failure  of  composite 
laminates  usually  occurs  before  catastrophic  disbond  propagation 
is  initiated  under  quasi-static  load,  although  slow  disbond 
growth  may  occur  under  cyclic  loading  (ref.  1).  Also,  when  the 
disbond  length  is  only  an  order  of  magnitude  higher  than  the 
diameter  of  fibers,  or  fiber  spacing,  use  of  elasticity  theory 
with  the  assumption  that  the  individual  layers  are  homogeneous 
becomes  questionable  and  further  refinements  in  analysis  seems 
either  impossible  or  too  complicated.  It  appears,  therefore, 
that  for  all  practical  purposes  shear  deformation  theory  used 
herein  is  quite  adequate  for  a  wide  range  of  disbond  size.  For 
cyclic  loading,  semi-empirical  growth  laws  of  the  type  given  be¬ 
low  are  commonly  used  for  metals,  where  a  single  mode  of  deforma¬ 
tion  (I,  II  or  III)  is  usually  dominant. 

C ( AK-AqK) n  (i) 

crack  length , 
number  of  cycles, 
empirical  constants 

stress  intensity  factor  range  =  K  -K  .  ,  and 

1  r  max  mm 

the  threshold  value  of  AK  below  which  propagation  does 
not  occur. 

Equation  (1)  is  commonly  used  for  studying  self-similar 
crack  growth  in  2-D  problems  and  has  been  used  for  modeling  the 
growth  of  disbonds  in  composite  laminates  (refs.  2,3).  Since 
the  model  used  here  does  not  yield  the  stress  intensity  factors, 
use  of  the  following  relationship  may  be  attempted. 


da  _ 
dN 

where  a  = 
N  = 
C,n  = 
AK  = 
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(2) 


a§  =  cl[&(/g)  -  A0(/^)]n 

where  A(/G)  =  /G„ ~ -  /G  !  ,  and 

max  min 

Aq(/G)  =  the  threshold  value  of  A ( /G)  below  which  propaga¬ 
tion  does  not  occur. 

da  can  be  considered  as  the  movement  of  the  disbond  front  in  the 
direction  of  current  normal.  Equation  (2)  is  consistent  with 
(1)  because  the  strain  energy  release  rate  is  proportional  to 
the  square  of  the  stress  intensity  factor. 

For  the  buckling  failure  outlined  in  case  (2) ,  use  of  an 
iterative  procedure  yields  critical  loads  based  on  linearized 
theory  of  elastic  instability.  Use  of  such  a  lower  order  theory 
for  sandwich  beams  containing  a  disbonded  compression  skin  has 
been  found  to  predict  buckling  failure  at  loads  close  to  exper¬ 
imentally  observed  values  as  long  as  the  disbond  length  is 
sufficiently  large.  For  small  disbond  lengths,  the  predictions 
tend  to  become  non-conservative  because  of  two  reasons;  namely: 
(i)  inelastic  buckling;  and  (ii)  inherent  limitations  of  shear 
deformation  beam  theory  for  small  length  to  depth  ratio. 
Theoretical  results  are,  however,  useful  for  a  wide  range  of 
disbond  length. 

METHOD  OF  ANALYSIS 

Double  Fourier  transform  is  made  use  of  (see  Appendix  A-2) 
to  handle  the  system  of  governing  differential  equations  given 
in  Appendix  A-l.  Three  unknown  functions  are  introduced  to 
characterize  the  discontinuities  in  displacements  in  three 
directions  over  the  area  of  the  disbond.  Elliptic  disbonds  are 
mapped  onto  a  circle  (r_<l  ,0<9j<2it)  by  a  simple  transformation 
and  the  unknown  functions  u*(r,0)  and  the  three  prescribed  trac¬ 
tions  t^(r,0),  (k=l,2,3)  are  expanded  in  Fourier  series  in  6. 
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u*(r,0)  =  Z  r|r|u*  (r")e' 

P  =  -°°  ** 


tk(r'0)  =  E  fckp(r)e 

P=-0O  r 


where  u*  ^  ^  and  t^(_pj  are  complex  conjugates  of  u£p  and  t^p , 

respectively • 

2 

Introduction  of  the  variables  u^  (r  )  and  t*  (r)  given  by 


r  v(r2’  *  ai  ukPlr2>  >  0<r<1 


Vr|  ‘  /  VrlrP+ldr 


and  expansion  of  real  and  imaginary  parts  of  u  and  t*  in  a 

2  ^P 

series  of  Legendre  polynomials  P^_^(2r  -1),  n=l,2...°°  as  well  as 
other  algebraic  manipulations  results  in  an  infinite  system  of 
linear  algebraic  equations  which  must  be  solved  to  determine  the 
unknown  constants  u  for  known  values  of  the  expansion  coef¬ 


ficients  t- 


kpq2n 


jmp^n'  . 


2  3  00  00  jmq-.n' 

Z  Z  Z  Z  Z,  iu,  =t.  j  =  1,2,3 

q.  =1  k=l  p=0  n=l  kpq2n  kpq2n  :mgin  _ 


The  expression  for  the  elements  of  the  matrix  l,  other  de¬ 
tails  and  method  of  calculation  of  strain  energy  release  rates 
Gj,  Gjj  and  G^^  are  given  in  Appendix  A-2.  For  the  purpose  of 
numerical  solution  it  is  necessary  to  truncate  the  system  (equa¬ 
tion  5)  such  that  m,p<m  and  n,n'<n  .  For  orthotropic  lam- 
inates  with  x,y  as  principal  axes,  the  set  of  equations  for 
ukpln  anc*  ukp2n  are  uncouPled  a^d  are  solved  separately. 

The  condition  that  the  system  of  equations  is  homogeneous, 
i.e.  the  determinant  of  the  coefficient  matrix  is  equal  to  zero, 
yields  the  critical  loads  for  the  buckling  problem.  Since  the 
values  of  prestress  enter  each  of  the  terms  of  the  coefficient 
matrix  in  a  complex  fashion,  the  condition  yields  a  complicated 
transcendental  equation  and  therefore,  the  critical  loads  must 
be  determined  by  trial  and  error  procedures. 
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EXPERIMENTAL  PROGRAM 


SHEAR  TESTS  ON  THICK  DISBONDED  LAMINATES 

Specimen  Geometry  and  Fabrication 

Thick  64-ply  Graphite/Epoxy  laminates  of  stacking  sequence 
of  [  (0 ^/±4  52/+45 2/0 s] s  were  tested  in  three  point  bending  with 
a  span  of  152.4  mm.  Specimens  were  254  mm.  in  length,  76.2  mm. 
in  width  and  approximately  9.73  mm.  in  thickness.  Two  teflon 
disbonds  were  implanted  at  the  midsurface  {between  plies  32  and 
33)  in  each  specimen.  The  two  disbonds  were  placed  symmetri¬ 
cally  with  respect  to  the  center  of  the  beam  span.  To  minimize 
the  influence  of  free  edges,  widths  of  the  flaws  did  not  exceed 
one-half  of  the  specimen  width.  Criticality  of  various  flaw 
geometries  were  examined  by  using  five  different  sizes  and 
shapes  of  disbond,  namely:  (i)  circular  disbonds  of  25.4  mm., 
31.75  mm.  as  well  as  38.1  mm.  in  diameter;  (ii)  oblong  shaped 
disbonds  with  semicircular  ends  of  dimensions  31.75  mm.  x  25.4 
mm.  and  38.1  mm.  x  25.4  mm.  (fig.  3).  The  oblong  shaped  dis¬ 
bonds  were  placed  in  the  beams  with  their  largest  dimension 
parallel  to  the  length  of  the  beam.  Centers  of  the  disbonds 
were  located  either  ±38.1  mm.  or  ±50.8  mm.  from  the  center  of 
the  beam.  The  Teflon-Teflon  disbonds  were  constructed  from  a 
stock  FEP  Teflon  film  .0254  mm.  in  thickness  by  doubling  the 
sheet  and  using  a  steel  punch  to  cut  it  in  the  shape  of  the  dis¬ 
bond  geometry  (fig.  4). 

Hercules  AS/3501-6  Graphite/ Epoxy  material  in  the  form  of 
304.8  mm.  wide  prepreg  tape  was  laminated  into  panels  measuring 
355.6  mm.  in  width  by  355.6  mm.  in  length.  Each  panel  was  as¬ 
sembled  from  two  half  panels  of  32  plies,  each  having  the  stack¬ 
ing  sequence  (04/±452/+452/04)  .  Four  specimens  were  obtained 

from  each  processed  panel.  After  being  processed,  the  panel  was 
inspected  ultrasonically  to  check  laminate  integrity  and  deter¬ 
mine  the  exact  location  of  the  specimens  within  the  processed 


panel.  Individual  specimens  were  cut  from  the  panel  using  a 
diamond  impregnated  saw. 

Test  Procedure 

The  thick  laminated  beams  were  tested  in  a  three-point  bend¬ 
ing  fixture,  shown  in  figure  5.  The  9.53  mm.  round  steel  load 
reaction  pins  served  to  transmit  a  line  load  across  the 

width  of  the  beam.  Support  arms  ©  (fig.  6)  were  rigidly  con¬ 
nected  to  the  main  frame  rail  ®  •  During  preliminary  tests 
the  9.53  mm.  loading  nose  ©  caused  localized  delaminations  and 
surface  damage  which  was  corrected  by  inserting  a  12.7  mm.  x  76.2 
mm.  x  3.18  mm.  copper  pad  between  the  steel  nose  and  specimen 

surface . 

Static  tests  were  conducted  in  an  Instron  model  TTCML  Univer¬ 
sal  testing  machine  using  a  crosshead  speed  of  0.05  cm. /min. 
Failure  was  determined  as  the  maximum  load  recorded  on  the  load 
deflection  curve.  Fatigue  loading  was  performed  in  an  Instron 
model  1321  servo-hydraulic  test  frame.  Sinusoidal  loading  was  ap¬ 
plied  at  a  frequency  of  10  Hz  with  a  minimum  to  maximum  stress 
amplitude  ratio  of  0.1.  The  test  duration  was  nominally  106  cy¬ 
cles  with  periodic  monitoring  of  the  flaw  propagation  at  selected 
intervals  during  the  load  history. 

Prior  to  testing,  specimens  were  ultrasonically  C-scanned  to 
provide  initial  flaw  measurements  and  locate  the  midpoint  between 
the  two  disbonds.  A  line  was  marked  on  each  beam  corresponding 
to  the  center  of  the  span.  The  specimen  was  fatigue  tested  for 
an  interval,  then  C-scanned.  In  the  beginning,  C-scans  were  taken 
after  every  5,000  or  10,000  cycles.  Once  a  feel  for  the  growth 
characteristics  was  obtained,  appropriate  test  intervals  were  used 
to  expedite  testing,  based  on  disbond  size  and  S-level.  Once 
propagation  of  the  defect  occurred,  intervals  were  reduced  for  de¬ 
tailed  characterization  of  flaw  growth.  Termination  of  a  test 
occurred  when  the  flaw  failed  to  propagate  despite  continued 


cycling,  or  when  the  specimen  failed  through  catastrophic  crack 
propagation.  Tests  were  conducted  at  S-levels  of  0.4,  0.5  and 
0.6  depending  on  the  size  of  the  disbonds. 

Pulse  echo  ultrasonic  inspection  was  used  to  nondestructive- 
ly  monitor  the  flaw  propagation  during  the  tests.  The  ultrasonic 
system  consists  of  a  model  5052  Panametrics  ultrasonic  analyzer, 
Tektran  tank  with  scanner  and  hot  pen  recorder,  and  a  Tektronix 
dual  trace  oscilloscope.  The  gated  peak  detector  section  of  the 
ultrasonic  analyzer  was  used  to  select  the  portion  of  the  wave¬ 
form  corresponding  to  the  implanted  defects  enabling  the  study 
to  focus  on  damage  propagation  from  the  implanted  flaws. 

Test  Results 

The  results  of  the  static  tests  were  primarily  utilized  in 
the  determination  of  the  stress  levels  used  in  the  fatigue  pro¬ 
gram.  The  data  are  presented  in  table  1.  All  specimens  first 
failed  in  flexure,  then  subsequently  failed  in  shear  (see  fig.  7) . 
Observation  of  specimen  failures  indicates  tensile  failure  of  ±45° 
plies  adjacent  to  the  outer  0°  plies  (plies  59-60) .  This  led  to 
subsequent  failure  of  the  0°  plies.  Additional  loading  of  the 
specimen  caused  the  specimen  to  fail  in  shear. 

Load  levels  for  fatigue  tests  were  determined  based  on  a 
static  failure  load  of  29,400N,  which  was  estimated  from  the  first 
few  of  the  static  tests.  The  data  obtained  from  the  fatigue  test¬ 
ing  program  were  compiled  and  produced  as  tables  and  graphs  of 
crack  length  versus  cycles.  These  results  can  be  found  in  refer¬ 
ence  19.  C-scans  of  some  of  the  specimens  are  shown  in  figures  8- 
18. 

Two  disbonds  in  each  specimen  gave  rise  to  four  crack  fronts. 
Disbond  I  was  the  upper  disbond  (nearest  the  l.D.  number),  II  the 
lower  disbond.  A  and  B  were  symbols  used  to  denote  the  upper  and 
lower  halves  of  disbonds  I  and  II,  respectively.  Likewise,  C  and 
D  denoted  the  upper  and  lower  halves  of  disbond  II  (see  fig.  19) . 
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These  four  disbond  halves  A,  B,  C  and  D  initiated  the  four  crack 
fronts.  Measurements  of  crack  lengths  were  taken  at  three  points 
along  each  crack  front  (locations  1,  2  and  3  in  the  figure). 

These  measurement  locations  were  independent  of  disbond  geometry 
or  load  level. 

The  number  of  cycles  to  which  the  specimens  were  loaded  is 
listed  in  table  2,  where  the  defect  size  in  inches  indicates  the 
largest  dimension  of  the  initial  flaw.  The  letters  C  and  0 
indicate  circular  or  oblong  shaped  flaws.  For  all  oblong  shaped 
flaws  the  minimum  dimension  was  25.4  mm.  The  location  indicates 
the  distance  of  the  center  of  the  flaw  from  the  center  of  the 
beam. 


SANDWICH  BEAM  TESTS 

Sample  Geometry  and  Fabrication 

The  specimen  geometry  consisted  of  a  22.0  in.  x  3.0  in.  x 
1.6  in.  (55.9  cm.  x  7.6  cm.  x  4.1  cm.)  sandwich  beam  with  a  12- 
ply  0°  laminate  on  the  tensile  side  and  a  [(0/±45)  ]  laminate  on 
the  compressive  side.  The  core  material  was  24  lb.  Hexcel  alumi¬ 
num  honeycomb  and  the  laminates  were  made  of  AS/3501-6  Graphite/ 
Epoxy.  The  Graphite/Epoxy  face  sheets  were  bonded  to  the  honey¬ 
comb  core  using  FM300M  sheet  adhesive. 

Before  processing  in  the  autoclave,  the  face  sheets  measured 
0.064  in.  (1.63  mm.)  and  the  honeycomb  core  measured  1.503  in. 
(38.2  mm.).  After  processing  the  finished  specimen  measured 
1.630  in.  (41.4  mm.) .  This  indicates  that  the  core  was  com¬ 
pressed  during  the  processing  and  bond  line  thickness  could  not 
be  measured  by  a  sum  of  the  parts  technique.  The  core  compres¬ 
sion  ranged  from  0.001  in.  (0.025  mm.)  to  0.003  in.  (0.076  mm.). 
Bond  line  thickness  was  approximately  0.001  in.  (0.025  mm.)  as 
measured  by  microscopy.  During  the  lay-up  of  the  [ (0/±45)  ] 
face  sheets,  two-layer  Teflon  defects  were  imbedded  between  the 


midplies,  or  between  the  -45°  plies,  three  plies  from  the  exposed 
face.  The  defects  were  of  either  circular  or  oblong  geometries 
and  were  made  by  punching  two-ply  0.001  in.  (0.025  mm.)  thick 
Teflon  film  against  an  abrasive  surface.  The  abrasive  surface 
causes  the  two  layers  of  Teflon  to  bond  together  and  form  a 
Tef lon-to-Tef Ion  defect.  The  study  required  the  fabrication  of 
18  sandwich  beams  with  the  defect  sizes  and  locations  given  in 
table  3.  In  addition  to  those  listed,  one  beam  was  fabricated 
with  1.5  in.  circular  defects  implanted  between  the  midplies  of 
both  the  12-ply  0°  and  the  [0/±452/0]s  laminates.  A  schematic  of 
the  aluminum  honeycomb,  Graphite/Epoxy  sandwich  beam  and  the  load 
application  points  is  shown  in  figure  20. 

Pretesting  Ultrasonic  C-Scan  Inspection 

Prior  to  testing,  all  samples  were  ultrasonically  inspected 
to  determine  defect  size  and  location.  Using  the  ultrasonic  C- 
scans  as  a  reference,  strain  gages  were  applied  to  both  the  ten¬ 
sile  and  compressive  faces  at  the  center  of  the  defect.  A  repre¬ 
sentative  C-scan  of  each  type  of  specimen  is  shown  in  figures  21a 
through  21f. 


Testing  Procedure 

The  experimental  test  setup  is  shown  schematically  in  figure 
22.  Load  and  strain  information  was  put  into  Vishay  model  2120 
bridge  amplifiers  and  the  conditioned  output  was  graphed  on  two 
X-Y  plotters.  Load  versus  compressive  strain  was  monitored  on 
plotter  No.  2.  Loading  was  discontinued  if  buckling  was  observed 
or  upon  laminate  failure.  If  buckling  occurred,  the  sandwich 
beam  was  removed  from  the  test  machine  and  an  ultrasonic  C-scan 
performed  on  the  sample  to  determine  the  presence  and  extent  of 
any  propagated  interlaminar  damage.  After  the  ultrasonic  nonde¬ 
structive  inspection,  the  specimen  was  tested  beyond  the  buckling 
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load  and  removed  from  the  test  machine  at  intervals  to  character¬ 
ize  the  extent  and  growth  of  damage.  Figure  23  shows  ultrasonic 
C-scans  of  a  near  surface  delamination  specimen  at  various  stages 
Figures  24  through  28  are  representative  photographs  of  failures 
of  the  specimens.  Photographs  of  some  more  specimens  can  be 
found  in  reference  20. 


Test  Results 

Three  beams  of  each  of  the  previously  mentioned  configura¬ 
tions  were  tested  and  the  results  are  summarized  in  table  4.  For 
the  case  of  near  surface  defects,  which  showed  some  indication  of 
buckling  prior  to  final  failure,  the  test  data  are  reported  in 
table  5 . 

Typical  test  results  for  some  of  the  different  configura¬ 
tions  tested  are  given  in  graphic  form  in  figures  29  through  32. 
These  results  show  graphs  of  load  and  stress  versus  compressive 
strain.  Other  typical  results  can  be  found  in  reference  20.  In 
all  cases  involving  midply  defects,  the  failures  were  instantan¬ 
eous  and  catastrophic.  For  the  near  surface  defect  tests,  test 
results  indicated  that  failure  appeared  to  have  started  as  buck¬ 
ling  above  the  defect  region  and  progressed  with  increasing 
damage  to  the  surface  plies.  Failure  continued  with  the  instan¬ 
taneous  and  catastrophic  failure  of  the  plies  below  the  implanted 
delamination.  However,  except  for  one  specimen,  all  of  the  sam¬ 
ples  tested  showed  a  sudden  change  in  any  reversal  of  strain  in 
the  compression  skin.  The  last  of  the  samples  with  a  1.5  in. 
near  surface  defect  {see  table  5) ,  which  appeared  to  have  a  de¬ 
formed  shape  prior  to  testing,  showed  an  increase  in  compressive 
strain  with  initial  loading  and  gradual  decrease  at  a  later  stage. 
This  behavior,  which  is  the  characteristic  of  Euler  type  buckling, 
is  shown  in  figure  32. 
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ULTRASONIC  WAVE  PROPAGATION  STUDIES 


Sample  Geometries 

The  specimens  were  prepared  from  32-ply  laminates  of  AS/ 
3501-6  Graphite/Epoxy.  The  two  laminate  configurations  were 
[±450]  and  [0,o]  and  the  laminates  were  machined  into  1  in.  x 
9  in.  tensile  coupons.  Resin  samples  were  prepared  from  the 
AS/3501-6  resin  system  at  the  University  of  Wyoming.  The  speci¬ 
mens  were  cured  at  177°C  (350°F)  and  were  then  cut  into  2.54  cm. 
x  22.86  cm.  (1  in.  x  9  in.)  test  coupons.  The  surfaces  had  to 
be  ground  in  order  to  remove  surface  defects  which  consisted  of 
multiple  shallow  depressions  a  square  inch  or  more  in  area  and 
a  few  thousandths  of  an  inch  deep.  The  specimens  were  then 
post-cured  for  eight  hours  at  218°C  (425°F)  and  the  surfaces 
were  polished.  Occasionally  there  are  dark  flecks  in  some  of 
the  specimens.  These  appear  to  be  impurities  or  contamination 
in  the  resin  as  received  from  Hercules,  Inc.  The  grinding  and 
polishing  produces  an  excellent  surface;  however,  examination 
of  the  specimens  in  a  polariscope  indicate  some  residual  stress¬ 
es  which  were  not  completely  annealed  out  during  post-cure.  An 
attempt  at  a  204°C  (400°F)  post-cure  for  10  hours  was  almost 
totally  ineffective  in  reducing  the  residual  stresses. 

Signal  Processing  Facility 

The  signal  processing  facility  developed  for  the  analysis 
of  ultrasonic  waveforms  is  shown  schematically  in  figure  33. 

The  system's  input/output  (I/O)  components  are  a  DEC  LA36  paper 
terminal,  a  Tektronix  CRT  graphic  display  terminal,  a  paper  tape 
reader/punch,  and  an  Omnitec  telephone  modem.  These  components 
are  all  directly  connected  to  either  the  parallel  or  serial  in¬ 
terface  ports  of  a  Digital  LSI  11  minicomputer  which  was  expand¬ 
ed  by  an  additional  16  dual  card  slots  in  a  Digital  expansion 
unit . 
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The  data  sampling  and  digitizing  component  in  the  system  is 
a  Nicolet  digital  oscilloscope.  The  oscilloscope  has  a  maximum 
sampling  rate  of  20  MHz  and  digitizes  the  input  to  eight  bits. 

The  oscilloscope  has  the  ability  to  store  waveforms  on  a  self- 
contained  floppy  disk  and  can  also  output  waveforms  to  a  flatbed 
plotter.  Since  the  maximum  frequency  which  can  be  digitized  on 
this  unit  without  aliasing  is  10  MHz,  a  10  MHz  low-pass  filter, 
whose  3-dB  point  is  at  10  MHz  and  which  is  65  dB  down  at  11  MHz, 
is  used  to  pre-filter  all  inputs.  A  peripheral  interface  was 
constructed  to  join  the  Nicolet  oscilloscope  to  the  LSI  11  com¬ 
puter.  The  interface  expands  16  input  bits  and  25  output  bits 
on  the  LSI  11  to  112  input  bits  and  112  output  bits.  The  oscill¬ 
oscope  requires  25  input  bits  and  20  output  bits  for  control  and 
data  input  and  output.  The  remaining  bits  of  the  interface  will 
be  used  for  system  expansion.  All  output  bits  are  tristate-con¬ 
trolled  from  the  front  panel  of  the  peripheral  interface.  This 
allows  data  and  control  lines  to  be  set  in  the  proper  state  prior 
to  data  transmissions  and  also  allows  manual  operation  of  the  os¬ 
cilloscope.  In  order  to  perform  the  required  data  transmissions 
and  manipulations,  a  control  program  named  NICDAT  was  written  in 
LSI  11  machine  language.  The  details  of  this  program  and  a  list¬ 
ing  of  support  programs  are  given  in  reference  21.  The  primary 
use  of  the  NICDAT  program  is  to  control  the  peripheral  interface 
and  to  transfer  data,  between  set  limits,  from  the  Nicolet  oscil¬ 
loscope  to  the  Burroughs  B7700  computer  where  the  data  can  be 
stored  and  analyzed. 

The  support  programs  allow  the  user  to  make  velocity  of 
sound  measurements,  attenuation  measurements,  interpolate  the 
waveform  using  frequency  domain  information  and  obtain  frequency 
spectra.  The  velocity  of  sound  is  obtained  from  measurements  of 
the  phase  of  the  front  surface  echo  and  back  surface  echo  indepen 
dently.  The  phase  of  the  two  echoes  can  be  determined  from  the 
Fourier  transform  of  the  individual  echoes.  This  phase  informa¬ 
tion  can  be  used  to  determine  the  velocity  of  sound  in  the  sample 


as  a  function  of  frequency,  and  the  theory  behind  the  measurement 
is  diagrammed  in  figure  34.  A  phase  measurement  is  made  of  the 


first  surface 


>q  and  at  the  second  surface  .  The  thickness  of 


the  sample  must  be  measured  at  the  point  on  the  sample  where  the 
phase  measurements  are  made.  With  this  information,  the  velocity 
of  sound  can  be  evaluated  using  equation  (6) . 


c  (u)  = 


R7  (oj)  -  <J)q  (o>) 


The  results  of  this  computation  are  shown  for  a  6.35  mm.  (1/4  in.) 
Graphite/Epoxy  sample  in  figure  35. 

Attenuation  of  sound  measurements  are  made  using  either  fre¬ 
quency  domain  information  or  amplitude  information.  Using  the 
Fourier  transform  of  the  front  surface  echo,  back  surface  echo 
and  second  multiple  echo,  the  attenuation  coefficient  a  can  be 
evaluated  as  a  function  of  frequency  w  using  equation  (7) . 


a  (w)  = 


R3(w) 
R2  (w) 


R2  (w) 


R1(w) ,  R2(w)  and  R3 (w)  are  the  magnitudes  of  the  Fourier  trans¬ 
forms  of  the  three  echoes  mentioned  and  d  is  the  thickness  of  the 
sample.  An  example  of  this  measurement  for  32-ply,  0°  Graphite/ 
Epoxy  is  shown  in  figure  36.  The  vertical  lines  in  the  waveform 
at  the  top  of  the  figure  indicate  the  boundaries  of  the  Fourier 
transform  of  each  echo. 

The  attenuation  of  sound  may  also  be  measured  using  the  am¬ 
plitudes  of  the  echoes;  however,  a  in  this  case  will  not  be  a 
function  of  frequency  oj.  For  this  measurement,  equation  (8)  is 
used . 


(8) 


a  =  ^  In  [<R3/R2)  -  (R2/R1) ] 

where:  a  is  the  linear  amplitude  attenuation  coefficient 

-1 

in  l.JTl  ; 

d  is  the  sample  thickness  in  mm; 

R-^  is  the  amplitude  of  the  front  surface  echo; 

R3  is  the  amplitude  of  the  back  surface  echo; 

R^  is  the  amplitude  of  the  second  multiple  echo. 

Due  to  the  inversion  which  occurs  at  the  back  surface  interface, 
r2  and  R^  are  measured  on  the  opposite  side  of  the  axis  from  R  . 
The  transducer  used  in  all  of  the  velocity  and  attenuation  meas¬ 
urements  is  a  5  MHz  type .  The  attenuation  measured  at  5  MHz  for 
a  0°  sample  was  0.053  mm.  1  using  the  frequency  technique,  and 
the  attenuation  was  0.061  mm.  1  using  the  amplitude  technique. 
There  is  a  13%  difference  between  these  measurements  and  this 
difference  is  typical  of  Graphite/Epoxy  composites.  The  dif¬ 
ference  decreases  to  approximately  5%  for  more  homogeneous  mate¬ 
rials  such  as  the  neat  resin  samples  tested. 

Data  Scattering  Test  Results 


In  order  to  determine  the  reproducability  of  the  ultrasonic 
velocity  and  attenuation  measurements,  eight  sets  of  data  were 
taken  at  several  locations  in  a  32-ply  ±45°  Graphite/Epoxy  test 
coupon.  Four  measurements  (FORSCATA-FORSCATD)  were  made  at  the 
center  of  the  sample  without  moving  the  transducer  and  four  meas¬ 
urements  (FORSCATE-FORSCATH)  were  made  around  a  1.27  mm.  (0.050 
in.)  diameter  circle  90°  apart.  The  results  of  the  velocity 
measurements  are  tabulated  at  1  MHz  intervals  in  reference  21, 
where  the  data  are  also  plotted  graphically.  A  5  MHz  nonfocused 
transducer  was  used  to  make  all  measurements.  The  transducer  was 
positioned  above  the  sample  for  a  main  bang  delay  of  137  ysec. 
(approximately  10.16  cm.).  Although  efforts  were  made  to  keep 
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the  transducer  stationary,  it  is  indicated  by  the  data  that  the 
sample  moved  after  the  first  set  of  data  was  taken.  The  remain¬ 
ing  data  for  FORSCATB  through  FORSCATH  indicate  a  very  low  degree 
of  scatter  at  the  5  MHz  frequency,  as  can  be  seen  from  the  se¬ 
lected  statistics  given  in  table  6 . 

Surface  Texture  and  Post-Curing  Test  Results 

Measurements  were  made  to  determine  what  effect  either 
polishing  the  surface  of  the  graphite  sample  or  post-curing  the 
sample  would  have  on  the  velocity  of  sound  and  attenuation  meas¬ 
urements.  Initial  measurements  were  made  on  the  top  half  and 
bottom  half  of  a  32-ply  ±45°  Graphite-Epoxy  sample.  These  meas¬ 
urements  are  reported  in  reference  21.  The  top  half  of  the 
sample  was  then  polished  on  both  faces  to  remove  the  rough  sur¬ 
face  which  is  left  when  the  peelply  is  removed  after  curing. 
Velocity  and  attenuation  measurements  were  then  made  over  the 
smooth  half  of  the  specimen.  Comparing  values  of  velocity  and 
attenuation  before  and  after  polishing  indicates  a  very  slight 
change  in  velocity  and  no  consistent  change  in  attenuation  over 
the  entire  frequency  range.  The  sample  was  then  post-cured  at 
40 0 °F  for  eight  hours  and  velocity  and  attenuation  measurements 
were  made  over  both  the  smooth  and  rough  surfaces.  The  results 
show  no  consistent  change  in  the  velocity  of  sound  after  post¬ 
curing  and  only  a  slight  increase  in  attenuation  over  the  entire 
frequency  range.  Since  these  results  indicated  no  general  change 
due  to  polishing  or  post-curing,  the  remaining  samples  in  the 
test  program  were  left  rough  and  were  not  post-cured. 

Moisture  Conditioning  of  Neat  Resin  and 
Graphite/Epoxy  Samples 

Neat  resin  samples  which  were  2.54  cm.  x  0.508  cm.  x  22.86 
cm.  (1  in.  x  0.2  in.  x  9  in.)  and  Graphite/Epoxy  samples  which 
were  either  [±45gJs  or  1 0 0 ] were  moisture  conditioned  by 


submersion  in  water  at  90°C  .  The  samples  were  weighed  before  and 
after  the  four  month  con  <  >ning  process  and  the  percentage  in¬ 
crease  in  weight  due  to  .  cure  absorption  was  14.91%  for  the 
resin,  2.73%  for  the  ±45°  samples  and  2.67%  for  the  0°  samples. 
Ultrasonic  velocity  and  attenuation  measurements  were  made  before 
and  after  conditioning.  The  results  for  a  frequency  of  5  MHz, 
which  are  summarized  in  tables  7  and  8,  indicate  no  significant 
change  in  the  velocity  of  sound  for  Graphite/Epoxy  as  well  as 
neat  resin  materials.  The  attenuation  measurements  were  made  us¬ 
ing  both  frequency  and  amplitude  information  and  the  results  in¬ 
dicate  significant  changes  in  the  attenuation  coefficients  of  all 
materials.  Results  in  table  8  show  that  it  can  be  concluded  with 
95%  confidence  that  moisture  content  alters  the  attenuation  coef¬ 
ficients.  Values  for  attenuation  and  velocity  at  1  MHz  intervals 
for  all  the  samples  before  and  after  moisture  conditioning  as 
well  as  sample  plots  of  these  data  can  be  found  in  reference  21. 

In  preparation  for  fatigue  testing,  static  tests  were  made 
on  the  ±45°  samples  and  the  neat  resin  samples.  Modulus  informa¬ 
tion  was  also  taken  on  the  neat  resin.  The  dry  resin  had  an 
average  ultimate  stress  of  30.63  MPa  (4442  psi) ,  an  average 
Poisson's  ratio  of  0.377,  and  an  average  Young's  modulus  of 
3.576  GPa  (5.187  x  105  psi).  The  dry  ±45°  graphite  samples  were 
tested  for  ultimate  stress  and  indicated  an  average  of 

172.9  MPa  (25,073  psi)  for  four  samples.  After  moisture  condi¬ 
tioning,  both  resin  and  ±45°  samples  were  tested  for  residual 
strength.  The  results  for  ultimate  stress  were  19.90  MPa  (2886 
psi)  for  the  resin  and  161.4  MPa  (23,415  psi)  for  the  ±45°  graph¬ 
ite  samples.  For  the  resin  samples,  all  failures  were  in  the 
center  of  the  test  section  and  figures  37  and  38  show  the  frac¬ 
tured  ends  of  the  dry  and  moisture  conditioned  resins. 

Moisture  conditioning  also  caused  a  significant  darkening  of 
the  color  of  the  resin  samples,  which  is  not  evident  in  the  photo¬ 
graphs.  Ultrasonic  C-scans  were  made  of  statically  tested  dry 
and  moisture  conditioned  ±45°  samples.  The  C-scans  indicate  a 


much  larger  degree  of  damage  outside  the  failure  region  for  the 
moisture  conditioned  samples.  These  C-scans  are  shown  in  fig¬ 
ures  39  and  40.  The  large  amount  of  damage  present  in  the 
moisture  conditioned  samples  was  not  present  prior  to  testing. 

Fatigue  Conditioning  of  [±45g]s  Graphite/Epoxy 

_ and  Neat  Resin  Samples _ 

Neat  resin  samples  were  fatigued  at  a  50%  S-level  for 
200,000  cycles  at  a  5  Hz  rate.  Velocity  and  attenuation  meas¬ 
urements  were  made  at  50,000  cycle  intervals  and  the  results 
indicate  no  changes  in  either  measurement. 

The  [± 4 5 _ ]  Graphite/Epoxy  samples  were  fatigued  at  75%, 
65%,  and  50%  S-levels.  The  samples  tested  at  75%  S-level  failed 
after  an  average  of  9030  cycles.  The  samples  cycled  at  S=65% 
failed  after  78,540  cycles  and  those  fatigued  at  S=50%  had  not 
failed  after  3,000,000  cycles.  Velocity  and  attenuation  meas¬ 
urements  were  made  on  the  samples  before  testing  and  after  fail¬ 
ure  as  well  as  at  10  intervals  during  testing.  Select  statis¬ 
tics  for  measurements  made  at  four  locations  for  all  load  cycle 
intervals  are  summarized  in  table  9.  Results  of  a  regression 
analysis  performed  for  examining  the  possibility  that  some  cor¬ 
relation  exists  between  the  measured  properties  and  number  of 
load  cycles  are  given  in  table  10.  It  appears  that  no  signifi¬ 
cant  correlation  can  be  established.  The  same  conclusion  is 
obtained  from  the  results  of  a  paired  "t"  test  reported  in 
table  11.  Ultrasonic  C-scans  were  performed  on  the  samples 
after  testing  or  failure  and  are  shown  in  figure  41.  The  scans 
indicate  no  damage  outside  the  failure  zones  and  no  damage  at 
all  after  3,080,000  cycles  at  50%  S. 
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ANALYTICAL/EXPERIMENTAL  DATA  CORRELATION 


DISBONDS  IN  SHEAR  ENVIRONMENT 


Strain  energy  release  rates  along  the  periphery  of  elliptic 

disbonds  of  various  sizes  located  at  the  midsurface  of  a 

[ (0^/±45^/+452/0^) g] s  graphite/epoxy  laminate  were  calculated  by 

the  use  of  the  methods  presented.  Elastic  properties  of  the 

unidirectionally  reinforced  composite  (0°  layer)  as  well  as 

stiffnesses  of  the  (0 ^/±45^/+452/0 laminate  are  listed  in 

table  12.  Since  D..  and  D„,  are  negligible  as  compared  to  D.. ..  , 

1  b  2.  o  11 

D^2  an<*  D22 '  bending-twisting  coupling  were  neglected  and  the 
two  (0<j/±452/+452/04 )  laminates  above  and  below  the  disbond 
were  considered  orthotropic.  Extensional  and  bending  stiffness¬ 
es  were  calculated  by  the  use  of  laminated  plate  theory  (equa¬ 
tions  A-1.4)  and  the  shear  stiffnesses  were  computed  with  the 
help  of  equations  (A-1.5)  and  (A-1.6)  given  in  Appendix  A-l. 

Shear  correction  factors  k^5  and  k^4  were  taken  equal  to  0.833. 

Retention  of  four  terms  (p,m  £  4  in  equation  A-2.57)  in  the 
Fourier  expansions  in  9  (equations  A-2.33,  A-2.34  and  3)  and 
four  terms  (n,n'  £  4  in  equation  A-2.57)  in  the  series  expansion 
in  r  were  found  to  yield  good  accuracy  in  obtaining  the  results. 
16  integration  points  were  sufficient  for  numerical  evaluation 
of  the  definite  integrals  with  respect  to  <j>  in  equations  (A-2.61) 
and  (A-2.62).  32  integration  points  were  found  to  be  enough  to 

yield  reasonable  accuracy  for  evaluation  of  the  infinite  inte¬ 
grals  with  respect  to  t,  in  equation  (A-2.62). 

Figures  42  and  43  show  the  variation  of  strain  energy  re¬ 
lease  rates  and  5j.-j.-j.  with  9  (along  the  disbond  periphery)  for 

various  disbond  geometries,  due  to  the  application  of  self-equil¬ 
ibrating  uniform  shear  stress  in  x-directions  on  the  disbond  sur¬ 
faces.  As  the  dimension  of  the  semi-axis  (in  x-direction)  is 
increased  keeping  the  other  semi-axis  L2  fixed,  values  of  Gjj 
first  increase  and  then  start  decreasing  whereas  values  of  G 
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go  on  increasing.  This  phenomenon  is  consistent  with  what  is 
observed  in  the  three-dimensional  elasticity  solution  of  an 
elliptical  crack  in  an  infinite  isotropic  medium  (ref.  22). 

Change  in  at  x=L^ ,  y=0  for  increasing  with  L2  held  con- 

stci.4t  is  more  clearly  demonstrated  in  figure  44.  Figure  45  shows 
the  variation  of  G^  at  x=L^ ,  y=0  with  L^/L^  as  L2  is  increased 
with  fixed.  It  should  be  noted  that  in  the  limiting  case  of 
Ll/L2~*0,  Pr°b-*-ern  reduces  to  that  of  generalized  plane 

strain.  Results  from  the  two-dimensional  elasticity  solution 
(ref.  2)  for  this  problem  are  also  shown  in  figure  45.  Although 
the  results  for  the  limiting  case  cannot  be  obtained  by  the 
methods  us^d  in  this  study,  the  results  in  figure  45  indicate 
that  laminated  plate  theory  solution  yields  values  of  /G 
which  are  about  7%  to  10%  lower  than  the  elasticity  solution. 
Similar  differences  were  observed  between  the  elasticity  solu¬ 
tion  (ref.  2)  and  the  strength  of  materials  type  solution 
(ref.  1)  as  indicated  in  figure  37  of  reference  2. 

Experimental  data  from  thick  beam  specimens  tested  under 
guasi-static  loading  do  not  indicate  any  failure  due  to  propa¬ 
gation  of  disbonds.  Therefore  it  was  not  possible  to  make  any 
assessment  of  critical  values  of  strain  energy  release  rates  to 
cause  quasi-static  propagation  of  disbonds.  Analytical  results 
indicate  that  for  a  disbond  of  38.1  mm.  diameter,  G ^  at  6=0° 
(x=19.05  mm.,  y=0)  can  reach  a  value  of  716  N/m  for  an  applied 
load  of  20,300  Newtons,  which  is  less  than  the  average  failure 
load  of  30,740  Newtons.  Tests  referred  to  in  reference  1  on 
25.4  mm.  wide  beams  containing  delaminations  of  the  same  width 
but  of  various  lengths,  indicated  a  critical  value  of  G^  of 
716  N/m  (see  ref.  2) .  In  the  specimens  used  in  the  current 
study,  values  of  G  ^  and  G  ^  at  points  other  than  at  9=0°  are 
significantly  lower  than  716  N/m  for  the  same  load  of  20,300 
Newtons.  In  fact,  G.^  is  a  maximum  at  0=0°  and  reduces  to  zero 
at  G=it/2  and  Gjjj  has  a  maximum  value  of  189  N/m  at  0=ir/2  and 
reduces  to  zero  at  9=0°  (see  figs.  42  and  43).  Therefore,  the 
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load  required  to  cause  disbond  propagation  is  likely  to  be  much 
higher  than  the  load  which  caused  flexural  failure  in  the  speci¬ 
mens.  As  reported  in  reference  19,  the  calculated  value  of 
load  to  cause  failure  of  the  45°  plies  is  32,300  Newtons,  which 
is  a  little  higher  than  the  average  experimental  value  of 
30,740  Newtons. 

Growths  of  disbonds  under  cyclic  loading  were  measured  and 
are  tabulated  in  reference  19.  The  results  show  that  the  growth 
in  y-direction  at  0=tt/2  is  not  at  all  significant.  Growth  rates 
in  x-direction  were  measured  at  six  points  (y=G  and  y=±9.53  mm.) . 
The  distance  of  the  two  disbond  fronts  at  y=0  from  the  center 
of  each  of  the  two  flaws  were  plotted  against  the  number  of 
cycles  N  for  each  specimen.  Representative  plots  are  shown  in 
figures  46-51.  Growth  rates  AL^/AN  were  then  computed  at  various 
points  along  the  L^  vs.  N  plots.  Values  of  at  these  points 
were  then  used  to  calculate  A  by  using  the  analytical  meth¬ 

ods  outlined  in  Appendix  A-2 ,  with  the  assumption  that  the  dis¬ 
bond  is  elliptic  with  semi-axes  L^  and  L2 .  Values  of  c^L^/dN  are 
plotted  against  A(/G~)  in  figure  52.  Experimental  results  in¬ 
dicate  that  the  growth  rates  are  extremely  small  for  S=0.4  and 
practically  no  growth  occurs  up  to  a  million  cycles  for  disbonds 
with  an  initial  value  of  2L1  £  31.75  mm.  Therefore  there  appears 
to  exist  a  threshold  value  A^  ( of  Af/fT^T)  below  which  propa¬ 
gation  of  disbonds  do  not  occur.  The  lowest  value  of  A ( /5  )  for 

lll/2 

which  some  growth  was  observed  is  of  the  order  of  9.4  (N/m) 

Least  square  fits  were  therefore  attempted  to  obtain  the  empiri¬ 
cal  constants  and  n  in  relationships  of  the  form: 


_  1/2 

with  a  series  of  values  of  varying  from  0  to  9  (N/m) 

The  parameters  for  two  curves  which  appear  to  give  the  closest 

fit  to  the  data  are  given  below  in  Newton,  meter  units. 


Aq(/giz),  (N/m) 1/2 

CL  (N 

,m 

units ) 

n 

7.0 

0.678 

X 

M 

O 

1 

VO 

2.979 

8.0 

2.565 

X 

10~9 

2.518 

The  curve  corresponding  to  the  first  set  of  parameters  is 
also  shown  in  figure  52.  The  curve  corresponding  to  the  second 
set  practically  coincides  with  the  first  in  the  range  of  A(/G~) 
shown  in  the  figure.  It  should  be  noted  that  the  growth  law 
given  by  equation  (9)  with  the  parameters  given  above  is  likely 
to  be  applicable  under  generalized  plane  strain  conditions. 
Curve  fits  attempted  in  reference  2,  which  is  appropriate  under 
generalized  plane  stress  conditions,  yielded  a  value  of 

_Q 

=  1.382  x  10  and  n  =  2.776  with  the  threshold  value 
V^>  equal  to  zero.  These  values  of  C.^  and  n  lie  between 
the  values  given  in  the  table  above.  For  plane  stress  condi¬ 
tions  effects  of  free  edges  are  likely  to  cause  the  flaws  to 
start  growing  under  much  smaller  values  of  A ( /G^^) .  Therefore, 
the  values  obtained  in  this  study  appear  to  be  consistent  with 
those  obtained  in  reference  2. 

BUCKLING  OF  DISBONDED  COMPRESSION  SKIN 


Properties  of  the  laminates  and  the  honeycomb  are  given  in 
table  13.  For  compression  skins  with  defects  in  the  midsurface 
there  is  no  extension-bending  coupling  and  twisting-bending 
coupling  is  not  too  large.  Therefore,  the  two  (0/+45/+45/0) 
laminates  can  be  considered  to  be  orthotropic.  For  compression 
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skins  with  defects  near  surface  (three  plies  below  exposed  sur¬ 
face)  there  exist  some  amounts  of  bending-extension  as  well  as 
bending-twisting  coupling.  For  the  present  study,  however,  all 
such  coupling  effects  were  neglected  to  reduce  the  computer  time 
needed  for  calculation  of  critical  loads  by  iterative  methods. 
Since  the  lowest  value  of  critical  load  was  to  be  determined, 
two  terms  in  the  Fourier  expansion  0  (p,m  £  2)  and  two  terms  in 
the  series  expansion  in  r  (n,n*  £  2)  were  retained  in  the  infin¬ 
ite  system  of  equations  (5).  12  integration  points  were  suffi¬ 

cient  for  numerical  evaluation  of  the  definite  integrals  with 
respect  to  <t>  in  equations  (A-2.61)  and  (A-2.62).  For  disbond 
sizes  not  too  large  compared  to  the  thickness  of  laminate  1,  32 
integration  points  were  found  to  yield  reasonable  accuracy  for 
evaluation  of  the  infinite  integrals  with  respect  to  c  in  equa¬ 
tion  (A-2.62).  For  larger  disbond  sizes  use  of  64  integration 
points  were  necessary. 

Figure  53  shows  the  variation  of  average  compressive  stress 
in  the  skin  required  to  cause  buckling  with  flaw  size  and  shape. 
The  pattern  is  similar  to  that  of  Euler  type  buckling.  The  cut¬ 
off  for  small  flaw  sizes  is  set  at  experimental  values  of  the 
average  compressive  stress  to  cause  failure  without  any  defect. 

Calculated  values  of  critical  stress  for  the  specimens  test¬ 
ed  are  compared  with  experimental  values  of  buckling  and/or 
failure  stress  in  table  14.  In  general,  the  experimental  values 
are  compar itively  much  higher  than  the  calculated  critical 
stress.  The  reason  for  this  discrepancy  appears  to  be  the  fact 
that  under  applied  load  the  sandwich  beams  have  certain  curvature 
which  increases  with  load,  but  the  portion  of  the  laminate  above 
the  disbond  cannot  have  a  buckled  shape  with  curvature  in  the 
same  direction,  because  of  constraints  imposed  on  it.  The  buckled 
shape  has  to  have  a  curvature  in  the  opposite  direction.  Studies 
on  buckling  of  initially  curved  bars  under  lateral  loads  (ref.  23) 
show  that  loads  to  cause  buckling  increase  rapidly  with  increas¬ 
ing  values  of  initial  deformation  and  sudden  reversal  of 
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deflection  occurs  at  the  critical  load.  Sudden  catastrophic 
failure  of  most  of  the  sandwich  beam  specimens  indicates  similar 
behavior.  It  may  be  noted  that  only  one  of  the  specimens  tested 
(specimen  No.  3  with  38.1  mm.  dia.  near  surface  defect)  showed 
a  gradual  change  in  the  deflection  pattern,  and  the  experimental 
value  of  average  compressive  stress  to  cause  strain  reversal  is 
lower  than  the  calculated  critical  stress  for  this  specimen.  A 
detailed  study  of  the  phenomenon  of  deflection  reversal  in  the 
types  of  specimens  tested  is  beyond  the  scope  of  the  present  work 
but  should  be  considered  in  future  studies,  since  the  simple 
Euler  type  buckling  assumption  underestimates  the  residual 
strength  of  sandwich  beams  tested  under  four-point  bending. 


DISCUSSIONS  AND  CONCLUSIONS 


In  the  present  study  a  NDE  methodology  for  assessing  critical¬ 
ity  of  disbond  type  defects  with  nearly  elliptic  two-dimensional 
planform  has  been  developed  and  its  usefulness  has  been  examined. 

In  addition,  effects  of  moisture  conditioning  and  cyclic  loading 
on  ultrasonic  wave  velocities  and  attenuation  have  been  studied. 

The  following  conclusions  and  comments  can  be  made  based  on  the 
results  obtained. 

1.  Linear  elastic  fracture  mechanics  approach  and  methods  of 
stress  analysis  based  on  a  modified  laminated  plate 
theory  yield  reasonably  accurate  values  of  strain  energy 
release  rates  along  the  periphery  of  a  disbond  in  a  lam¬ 
inated  composite  subjected  to  transverse  shear.  Based  on 
these  calculated  values  of  strain  energy  release  rates 
and  semi-empirical  crack  growth  laws,  growth  of  such  dis¬ 
bonds  under  cyclic  shear  can  be  predicted.  A  growth  law 
for  dominant  mode  II  behavior  under  generalized  plane 
strain  conditions  has  been  obtained  which  appears  to  be 
consistent  with  such  a  law  determined  in  reference  2  for 
generalized  plane  stress  conditions.  Further  studies  are 
needed  for  obtaining  growth  laws  under  combined  actions 
of  mode  II  and  mode  III  behavior.  Experimental  results 
indicate  that  for  the  types  of  specimens  tested  mode  II 
behavior  dominates  the  growth  of  the  disbonds.  By  the 
use  of  these  laws  and  appropriate  failure  models,  an 
assessment  of  residual  strength  and/or  residual  lifetime 
can  be  made . 

2.  Presence  of  disbonds  can  significantly  lower  the  load 
carrying  capacity  of  a  laminated  composite  under  compres¬ 
sion.  A  linearized  elastic  stability  analysis  can  be 
used  to  assess  criticality  of  disbonds  with  nearly  ellip¬ 
tic  planforms ,  but  the  analytical  predictions  appear  to 
underestimate  the  failure  loads  of  sandwich  beams  with 


delaminated  compression  skins  tested  in  four-point  bend¬ 
ing,  where  reversal  of  deflection  has  to  occur  before 
buckling  and/or  failure.  These  results  and  catastrophic 
failures  of  a  majority  of  sandwich  beams  tested  suggest 
the  need  for  further  studies  on  effects  of  deflection 
reversal  on  such  failures.  The  possibility  of  alterna¬ 
tive  test  plans  to  eliminate  such  effects  should  also  be 
examined  in  future  studies. 

3.  Moisture  conditioning  of  Graphite/Epoxy  composites  has 
some  effect  on  the  attenuation  of  sound  in  the  material, 
but  appears  to  have  little  or  no  effect  on  the  velocity 
of  sound.  Damages  induced  by  moisture  conditioning  are 
distributed  throughout  the  specimens.  Damages  caused  by 
cyclic  loading,  however,  are  highly  localized  near  the 
failure  sites.  Cyclic  loading  appears  to  cause  no  sig¬ 
nificant  changes  in  velocity  of  sound  and  attenuation 
prior  to  as  well  as  after  failure  of  Graphite/Epoxy 
specimens.  The  same  conclusion  holds  for  neat  resin 
samples  tested. 
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Table  1.  Laminated  Beam  Static  Test  Results 
([t,4/i452/i452/04]s)s 

Width=76.2  mm.,  Thickness=9 . 7  mm. 


Specimen 

Number 

Defect  Size 
(mm. /in . ) 

Failure  Load 
(N/lbs . ) 

1-1 

25.4/1.0 

25,310/5,690 

1-2 

25.4/1.0 

29,360/6,600 

1-11 

25.4/1.0 

36,740/8,260 

1-13 

25.4/1.0 

30 ,380/6 ,830 

1.25-1 

31.8/1.25 

34,790/7,820 

1.25-2 

31.8/1.25 

32,300/7,260 

1.5-2 

38.1/1.5 

30,380/6,830 

1.5-4* 

38.1/1.5 

26,690/6,000 

Average 


30,740/6,920 


Note:  Failed  in  Interlaminar  Shear 


\ 


Table  2.  Thick  Beam  Specimens  Tested  in  Fatigue 


Specimen  No.  Defect  S-Level 

t [ 0 4/  45-/  452/0.]  J  Size/Shape/Location 

s  (in.)  (in.) 


Total  Cycles 
Tested 


1-17 

1.0/C/1.5 

.  4 

<  625,000 

1-18 

1.0/C/l. 5 

.  4 

<1,075,000 

1.25-5 

1.25/0/1.5 

.4 

<1,000,000 

1.25/0/1.5 

.4 

<  725,000 

1.5-5 

1. 5/0/1. 5 

.4 

135,000 

1.5-6 

1. 5/0/1. 5 

.4 

115,000 

1-14 

1.0/C/l. 5 

.5 

180,000 

1-15 

1.0/C/l. 5 

.5 

130 ,000 

1.  25-3 

1.25/0/1.5 

.  5 

50,000 

1.25-4 

1.25/0/1.5 

.5 

100,000 

1.5-1 

1. 5/0/1. 5 

.  5 

27, 500 

1.5-3 

1. 5/0/1. 5 

.5 

25,000 

1-19 

1.0/C/l. 5 

.6 

32,000 

1-20 

1.0/C/l. 5 

.  6 

12,170 

1.25-7 

1.25/0/1.5 

.6 

11,500 

1.25-8 

1.25/0/1.5 

.6 

11,000 

1.5-7 

1 . 5/0/1 . 5 

.  6 

13,000 

1.5-8 

1. 5/0/1. 5 

.6 

16,000 

1.25C-1 

1. 25/C/1.5 

.5 

60,000 

1.25C-2 

1 . 25/C/l . 5 

.  5 

33,500 

1.5C-1 

1.5/C/1.5 

.5 

10,000 

1. 5C-2 

1.5/C/1.5 

.  5 

22,000 

1.25C-3 

1. 25/C/l. 5 

.6 

6,870 

1.25C-4 

L. 25/C/l. 5 

.6 

15,000 

/continued  next  page 
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Table  2  (continued) .  Thick  Beam  Specimens  Tested  in  Fatigue 


Specimen  No. 
[[04/±4S2/±452/04]s]s 


Defect  S- Level 

Size/Shape/Location 
(in. )  (in. ) 


Total  Cycles 
Tested 


1.5C-3 

1.5/C/1.5 

.6 

8,880 

1.5C-4 

1.5/C/1.5 

.6 

4,150 

1.0C2-1 

1.0/C/2.0 

.  5 

700,000 

1.0C2-2 

1.0/C/2.0 

.5 

328, 700 

1. 25C2-1 

1 . 25/C/2 . 0 

.  5 

118,000 

1.25C2-2 

1 . 2  5/C/2 . 0 

.5 

56,000 

•  1.5C2-1 

1.5/C/2.0 

.5 

19 ,840 

1.5C2-2 

1.5/C/2.0 

.5 

6,750 

1.0C2-3 

1. 0/0/2. 0 

.6 

65,650 

1.0C2-4 

1. 0/0/2. 0 

.6 

117,910 

1.25C2-3 
1. 25C2-4 


1.25/C/2.0 
1 . 25/C/2 . 0 


.6 

.6 


12,580 

47,090 


Table  3.  Specimen  Fabrication  Requirements  for 
Sandwich  Beam  Defect  Criticality  Study 


Defect  Size  Defect  Location  Number  of  Samples 


No  defect 

— 

3 

1.0-in.  Circular 

Near  surface* 

3 

1.5-in.  Circular 

Near  surface* 

3 

1.5-in.  Circular 

Midply 

3 

2.0-in.  Circular 

Midply 

3 

2.0-in.  2.5-in.  Oblong 

Midply 

3 

1.5-in.  Circular 

Midply-both  faces 

1 

TOTAL:  19 


*Near  surface  defects  are  located  three  plies  below  the 
exposed  surface. 


Table  4.  Sandwich  Beam  Failure  Test  Data 


Defect  Size 

Defect 

Location 

Maximum 

Load 

(kg) 

Maximum 

^T 

10  6 

Maximum 

Cc  6 
10”6 

No  defect 

2925 

5100 

15125 

No  defect 

— 

3425 

5775 

18100 

No  defect 

— 

3300 

5850 

15200 

1  in.  cir. 

Near 

2950 

— 

1  in.  cir. 

Near 

2750 

— 

+ 

1  in.  cir. 

Near 

2800 

— 

0) 
f— 1 

1.5  in.  cir. 

Near 

2900 

— 

o 

1.5  in.  cir. 

Near 

2730 

— 

g 

o 

1.5  in.  cir. 

Near 

2810 

— 

+ 

1.5  in.  cir. 

Mid 

2210 

3750 

9150 

1.5  in.  cir. 

Mid 

2715 

4325 

11300 

1.5  in.  cir. 

Mid 

1325 

2150 

5425 

2  in.  cir. 

Mid 

2275 

3925 

9850 

2  in.  cir. 

Mid 

2050 

3450 

8700 

2  in.  cir. 

Mid 

2500 

4225 

10500 

2  in.x2.5  in. 

Mid 

2350 

3800 

11725 

2  in.x2.5  in. 

Mid 

2435 

4025 

11000 

2  in.x2.5  in. 

Mid 

1960 

3450 

8425 

Table  5.  Sandwich  Beam  Buckling  Test  Data 


Defect  Size 

1  in.  cir. 

1  in.  cir. 

1  in.  cir. 

1.5  in.  cir. 

1.5  in.  cir. 

1.5  in.  cir. 


* 

Appeared  to 


Defect 

Location 

Buckling 

Load 

(kg) 

Buckling 

ET  6 

10"6 

Near 

1695 

2820 

Near 

1280 

2350 

Near 

140  J 

1770 

Near 

635 

1075 

Near 

775 

1330 

Near 

200* 

340 

have  a  deformed  shape 


Buckling 


6550 

5200 

4100 

2350 

2900 

105 


-38- 


Table  6.  Select  Statistics  for  Ultrasonic 
Property  Measurement 


5  MHz 


Standard  Average 

Location  Mean  Deviation  Variance  Deviation 

x  o  o2  £ | x-x | / n 


Velocity  (m. 

/sec . ) 

BCD 

3027.72 

0.4359 

0.1267 

0.720 

EFGH 

3028.27 

2.053 

4.2158 

1.635 

Attenuation 

(mm.  1) 

BCD 

0.093b 

0.0043 

1.25x10'  ' 

0.0022 

Ll’Gll 

0.0933 

0.0009 

3 . 58xl0~  5 

0.0049 

Table  7.  Select  Statistics  for  Moisture 
Conditioning  Results 


Mean 

x 


S  tandard 

Deviation  Variance 


Average 
Deviation 
E | x-x ] /n 


±45°  Graphite 


Velocity 

Dry 

(m./sec . ) 
3292.7 

54.13 

2511.9 

37.01 

Moist 

3192. 9 

80.52 

5557.3 

62.13 

Amp 

Atten 

Dry 

(mm.  ''') 
.06139 

. 01142 

.00011 

.00898 

Moist 

.09503 

.01076 

. 00010 

. 00817 

Freq 

Atten 

Dry 

(mm. 

.06897 

. 00949 

.00008 

.00682 

Moist 

.  11998 

.02670 

.00061 

.02118 

O'  Ciruphite 

Velocity 

Dry 

(m./sec . ) 
3213.4 

106.8 

9120 . 2 

77.68 

Moist 

3244.2 

87.1 

6065.8 

70.24 

Amp 

Atten 

Dry 

(mm.  ■*") 
.06244 

.00859 

.00006 

.00751 

Moist 

. 09812 

. 00749 

.00004 

.00586 

Fi'cq 

Atten 

Dry 

(mm .  ^ ) 
.07412 

.01841 

.  00027 

.01602 

Moist 

.12262 

.01945 

.00030 

.01366 

Neat  Resin 

Velocity 

Dry 

(m./sec. ) 
2921.7 

50.  36 

2173.6 

35.13 

Moist 

2869 . 4 

35.93 

1106.5 

28.60 

Amp 

Atten 

Dry 

(mm. 

. 12844 

.00529 

. 00002 

.00397 

Moist 

.17773 

.00337 

.00001 

.00258 

Freq 

Atten 

Dry 

(mm. 

. 12341 

.00774 

.00005 

.00567 

Moist 

. 18567 

.01585 

.00021 

.01420 

•nu  pi  mum 


-<|  V  V 
■c  >.  V  o 
oJ  ^  < 


u  O  •  • 

O  o  -•  O’ 
-i  a  ji 
•  ai  E  w 

o  >  <  u. 


4  ^  Q.  01 

«  a»  e  u 

z  >  <£  U. 


Table 

9.  Select 

Statistics  for 

±45°  Fatigue 

Results 

Position 

Mean 

X 

Standard 

Deviation 

0  i 

n-1 

Variance 

2 

a 

Average 

Deviation 

Z | x-x | /n 

Velocity 

f  (m./sec.) 

1 

3202.8 

82.45 

6042.5 

63.26 

4 

3234.2 

116.58 

12231.8 

89.44 

5 

3318.1 

120.29 

13022.3 

82.14 

6 

3286.2 

95.88 

8273.8 

69.40 

Attenuation  (mnu 


1 

.07366 

.OHIO 

.00011 

.00840 

4 

.05242 

.03633 

.00117 

.02425 

5 

.09871 

.01593 

.00023 

.01804 

6 

.07870 

.03835 

.00132 

.02512 

2- 


E(xjC)  -  b  +  mC ,  C  in  10  Cycles 


Table  11.  Paired  "t"  Test  for  t45°  Fatigue  Results 


Thousands  of  Cycles 


Table  12.  Elastic  Properties  of  0°  Layer  and 

Stiffnesses  of  (0 ./±45„/+45_/0  ) 

_  .  ,  4  2  2  4  s 

Laminate 

Transversely  Isotropic  0°  Layer 

Ea  =  125  GPa 

E  =  10  GPa 
T 

vA  =  0.28 

G„  =  5.8  GPa 
A 

(04/±452/+452/04) s  Laminate 

A-^  =  4.06  x  108  N/m. 

A12  =  0.79  x  108  N/m. 

A22  =  1.24  x  108  N/m. 

A,,  =  0.95  x  108  N/m. 
b  b 

A16  =  A26  =  °*° 

Ba£  =  °-°  (“'8 

K,.  _  =  2.12  x  107  N/m. 
bb 


D11  - 

877 

Nm 

°12  = 

131 

Nm 

°22  = 

218 

Nm 

°66  = 

160 

Nm 

°16  = 

°26 

= 

=  1,2,6) 

K44  =  1.75  x  107  N/m 


Ply  Thickness  =  0.1519  mm. 


Table  13.  Stiffnesses  of  the  Laminates  Comprising 
the  Compression  Skin  and  Properties  of 
the  Core  in  Sandwich  Beams 


(0/±45/+45/0)  Laminate,  Ply  Thickness  =  0.1354  mm. 


A11 

= 

5.64 

X 

107 

N/m. 

B11 

= 

0 

D11 

= 

4.51 

N.m. 

A12 

= 

1.68 

X 

107 

N/m. 

S12 

= 

0 

°12 

= 

0.48 

N.m. 

A22 

= 

2.50 

X 

107 

N/m. 

B22 

= 

0 

D22 

= 

0.86 

N.m. 

A66 

= 

1.92 

X 

107 

N/m. 

B66 

= 

0 

D66 

= 

0.62 

N.m. 

A16 

= 

0 

B16 

= 

0 

°16 

= 

0.29 

N.m. 

A26 

= 

0 

B26 

= 

0 

°26 

= 

0.29 

N  .m. 

K55 

= 

0.34 

X 

107 

N/m. 

K45 

= 

0 

K44 

= 

0.30 

x  10 

(0/±45)  Laminate,  Ply  Thickness  =  0.1354  mm. 


11 

= 

2.82 

X 

107 

N/m. 

B11 

= 

-1.55  x  103  N 

D11 

= 

0.46  N.m. 

12 

= 

0.84 

X 

107 

N/m. 

B12 

= 

0.49  x  103  N 

°12 

= 

0.093  N.m. 

’22 

= 

1.25 

X 

107 

N/m. 

B22 

= 

0.57  x  10 3  N 

°22 

= 

0.15  N.m. 

66 

= 

0.96 

X 

107 

N/m. 

B66 

= 

0.98  x  103  N 

°66 

= 

0.11  N.m. 

— 

0 

B_  c 

_ 

-0.53  x  103  N 

D->  c 

= 

-0.072  N.m 

16 

16 

16 

= 

0 

B.,  - 

— 

-0.53  x  103  N 

D„ 

_ 

-0.072  N.m 

26 

26 

26 

55 

= 

0.17 

X 

107 

N/m. 

K45 

= 

0 

K.  . 
44 

= 

0.15  x  107 

N/m. 
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Table  13  (continued) .  Stiffnesses  of  the  Laminates 

Comprising  the  Compression  Skin 
and  Properties  of  the  Core  in 
Sandwich  Beams 


(+45/0/±45/+45/0)  Laminate,  Ply  Thickness  =  0.1354  mm. 


A11  = 

8  , 

.46 

X 

107 

N/m . 

A12 

2. 

.52 

X 

io7 

N/m. 

A22 

3. 

.75 

X 

io7 

N/m. 

> 

03 

03 

II 

2. 

.88 

X 

io7 

N/m. 

A16  " 

0 

A26 

0 

K55 

0. 

,51 

X 

io7 

N/m. 

BX1  =  1.55  x  103  N 
B  =  -0.49  x  103  N 
B22  =  -0.57  x  103  N 
Bg6  =  -0.98  x  103  N 
B  6  =  0.53  x  103  N 
B26  =  0.53  x  103  N 

K.,  =  0 
45 


D  =  10.69  N.m. 

D12  =  3.05  N.m. 

E>22  =  4.58  N.m. 

D,,  =  3.50  N.m. 

66 

D  =  -0.22  N.m. 

D26  =  -0.22  N.m. 

K44  =  0.45  x  io'  N/m 


Core  Properties 


Thickness  =  38.17  mm. 

Young's  Modulus  =  5.17  GPa 
Shear  Modulus  (L)  =  2.55  GPa 
Shear  Modulus  (W)  =  0.62  GPa 
Subgrade  Modulus  in  x-Direction 
Subgrade  Modulus  in  y-Direction 
Subgrade  Modulus  in  z-Direction 


66.8  GN/m3 
16.3  GN/m3 
135.5  GN/m3 
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.001 '  THICK  TEFLON  FILM 


Figure  4.  Disbond  Fabrication  Method 


Test  Fixture 


Static  Interlaminar  Shear  Failure  (End  View) 


SPECIMEN  I.0C2-2 


Ficjure  8.  C-Scan  of  Specimen  1.0C2-2 


1 


SPECIMEN  I.5C-I 


SPECIMEN 


SPECIMEN  1.25  C-4 


SPECIMEN  15C-4 


Figure  18.  C-Scan  of  Specimen  1 . 5C- 
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TEFLON 
DISBOND  "I" 


MEASUREMENT 
LOCATION  LINES 


TEFLON 
DISBOND  "H" 


SPECIMEN 
D.  NO. 

CRACK  FRONT  "A" 

CRACK  FRONT  "B" 


iuii 


CRACK  FRONT  "C 


CRACK  FRONT  "D" 


DETECTED 
FLAW  GROWTH 


Figure  19.  Measurement  Locations  for  Four  Crack  Fronts 
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Figure  21f.  2"  x  2.:"  ■•  •;  Midply  Defect 


Figure  21  (continueo) .  C-I^iu; 


mdwich  Beam  Specimens 


i 


t 


I 


Instron 

Test  Machine 


Load  Data 


Strain 

Data 


Bridge 
1  Amplifier 


Load  Data _ 

_ «T 


r . . 

=  l 

si 

■ 

1 1 

t  k 

«c 

X-y  PLOTTER  I 
Load  vs.  compressive  strain 


X-y  PLOTTER  2 
Tensile  vs.  compressive  strain 


Figure  22.  Sandwich  Beam  Testing  Schematic 


1.5"  Circular  Near  Surface  Delamination  Prior  to  Testing 


1.5"  Circular  Near  Surface  Delamination  After  Buckling 


1.5" 


Circular  Near  Surface  Delamination 
After  2000  Kg  Applied  Load 


Figure  23.  Near  Surface  Delamination  Propagation  Sequence 
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Figure  26.  Failure  of  Sandwich  Beam  Composite 
with  a  1.5"  Circular  Defect 
Located  at  the  Near  Surface 


Figure  27.  Failure  of  Sandwich  Beam  Composite 
with  a  2.0"  Circular  Defect 
Located  at  the  Midply 


Figure  28.  Failure  of  Sandwich  Beam  Composite 
with  a  2.0"  Defect  Located  at  the 
Midply 
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Compressive  Strain  in  Sandwich  Beam  with  No  Defect,  Sample 


STRESS  VS.  €c 
LOAD  VS.  €C 


(6X)  OW)l 


80 


Figure  31.  Load  vs.  Compressive  Strain  in  Sandwich  Beam 
with  1"  Defect  Near  Surface,  Sample  1 


Defect  Near  Surface,  Sample  3  Crosshead  Speed=.2  cm. /min 


Figure  32.  Load  vs.  Compressive  Strain  in  Sandwich  Beam 
with  1.5"  Defect  Near  Surface,  Sample  3 


SIGNAL  PROCESSING  FACILITY 
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Figure  33.  Signal  Processing  Facility 
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P  ~  PATH  LENGTH  OF  SOUND 
X  =  WAVE  LENGTH 
<£0=  PHASE  AT  INTERFACE  0 
<£,  *  PHASE  AT  INTERFACE  1 
k  =  NUMBER  OF  CYCLES 
C  =  VELOCITY  OF  SOUND 

d  *  THICKNESS  OF  SAMPLE 

f  a  FREQUENCY  OF  SOUND 

W  «  FREQUENCY  IN  RAD/SEC 

x  =  c/f  c  =  \  f  w  5  2  ir  f 

--  Zwk  kX  =  p  p-. 

-g-f  2  y  .  P  u; 

c  tw }  =  UltL- 


34.  Phase  Velocity  of  Sound 


Figure  39.  C-Scans  of  Statically  Tested  Dry  and  Moisture  Conditioned 
32-Ply  ±45°  Graphite/Epoxy  Samples 


Figure  40.  Ultrasonic  C-Scans  of  Three  Moisture  Conditioned 
and  Statically  Tested  [+-45  ]  ,  Graphite-Epoxy 
Samples  s 


Figure  43.  Variation  of  GIIT  with  Flaw  Dimensions  for  [ (04/±452/+452/04) g] g  Laminate 


1  1.2  1.4  1.6  1.8  2  2.2  2.4 


i 


Figure  44.  Variation  of  GZI  at  (x=L1,y=0)  with  L^/I^ 
for  Fixed  Values  of  L2  for 
[ ( 0^/ ±4 5 2/+4 52/0 ^ ) s ^  s  Laminate 


for  Samples  1.5-7  and 


Figure  47.  L.  vs.  N  for  Samples  1.25-7  and  1.25- 


Figure  49.  L,  vs.  N  for  Samples  1.5C-1  and  1.5C- 


25- 


Figure  50.  L.  vs.  N  for  Samples  1.25-3  and  1.25- 


./cycle 


"P 


(0/+45/+45/0) 
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APPENDIX  A-l 


LAMINATED  PLATE  THEORY  WITH  EFFECTS  OF  SHEAR 


DEFORMATION,  SURFACE  TRACTIONS  AND  PRESTRESS 


The  following  displacement  field  is  usually  assumed 


in  a  laminated  plate  [ref.  7]: 


ux  =  u(x,y)  +  ztj>x(x,y) 
uy  =  v(x,y)  +  z<|jy(x,y) 


(A-l.l) 


u_  =  w(x,y) 


The  equilibrium  equations  in  terms  of  the  stress 

resultants  N _ ,  N  ,  N  ,  M  ,M  ,M  ,0,0  are  as 

xx'  yy'  xy'  xx'  yy'  xy'  wx'  uy 

follows : 


Nxx,x  +  Nxy,y  +  fcx  +  fcx  ~  0 


N  +  N  +  t  +  t  =  0 

xy,x  yy ,y  y  y 


Mxx,x  +  Mxy,y  +  I(tx  "  fcx*  °x  0 

M*y.*  +  Myy,y  +  f(ty  ‘  V  '  Qy  ‘  0 

Qx,x  +  Qy,y  +  Nxx  w,xx  +  2Nxy  w,xy 

+  N°  W  +  t+  +  t  =0 

yy  ,yy  z  z 

±  +  ± 

where  tx,  ty,  tz  are  the  surface  tractions  in  positive  x,  y 
and  z  directions  acting  on  the  surfaces  z  =  +|  of  the  plate, 
as  shown  in  figure  2.  To  study  the  possibility  of  buckling 
under  the  influence  of  the  in-plane  stress  resultants,  the 


(A-l. 2) 


r 


problem  is  considered  as  that  of  incremental  deformations. 

Thus,  all  quantities  of  second  order  are  omitted.  N°x, 

N°  and  N°  are  the  values  of  the  prestresses,  which  are 
xy  yy  c 

independent  of  spatial  coordinates  x  and  y.  Their  critical 
values  are  to  be  determined  such  that  an  alternative 
deformation  configuration  different  from  the  present  one 
is  possible. 

Substitution  of  the  relationships  between  the  stress 
resultants  and  derivatives  of  the  five  displacement  compo¬ 
nents  u,  v,  iji  ,  <p  and  w  results  in  the  following  system 
x  y 

of  differential  equations: 


AllU,xx  +  A66u,yy  +  2A16u,xy  +  A16v,xx 


+  A26V,yy  +  (A12  +  A66)v,xy  +  Bll’JJx,xx 


+  B,^\b  +  2B,  +  B. 

66vx,yy  16^x,xy  16^y,xx 


+  B-,.4 i  +  (B.  -  +  BcC)  ip 

26Yy,yy  12  66'ry,xy 


+  *x  +  fcx  =  0 


A. ,U  +  A~,U  +  (A..  +  A^ , ) U 

16  ,xx  26  ,yy  12  66  ,xy 


(A-1.3) 


+  A,  ,v  +  A_ »v  +  2Aocv  +  B.  ,\p 

66  ,xx  22  ,yy  26  ,xy  16vx,xx 


+  B~,i!/  +  (B,  _  +  B£C)iI)  +  B,ri)i 

26yx,yy  '  12  66,yx,xy  66vy,xx 


+  B22^y  ,yy  +  2B26*y,xy  +  fcy  +  t"  =  0 


$  -103- 
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B11U,XX  +  B66u ,yy  +  2B16u,xy  +  B16v,xx 
+  B26v ,yy  +  (B12  +  B66)v,xy  +  Dll^x,xx 
+  D66^x,yy  +  2D16^x,xy  +  D16^x,xx 
+  D26*y,yy  +  (D12  +  D66^yfxy  “  K55(*x  +  w,x) 
'  K45(*y  +  “,y>  ,+  T(tx  ‘  Si>  '  0 


B16u,xx  +  B26u ,yy  +  (B12  +  B66)u,xy 


+  B66Vfxx  +  B22V,yy  +  2B26V,xy  +  °16^x,xx 


+  D26^x,yy  +  (°12  +  D66)*x,xy  +  D66*yrxx 


+  D22^y,yy  +  2D26^y,xy  "  K45 (*x  +  w,x) 
-  K44»y  +  “,y>  *  f«=y  -  V  ’  ° 


IK55  +  N°x>",xx  +  2(K45  +  Nxy’",xy 


*  K44W,yy  *  K55*x,x  +  K45(*x,y  +  *y.x> 


+  K .  A\ii  +  t  +  t  =  0 

44yy,y  z  z 


The  laminate  stiffnesses  Bag '  Dag  = 

and  K  D  (a, 3  =  4,5)  are  calculated  from  the  plane  stress 
exp 

coefficients  Q  D  of  all  the  layers  in  the  plate  by  the 

ap 

following  relations: 
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H/2 


ag '  Bag  '  Dag 


I 


-H/2 


Qct3  {1'z'z‘ 


)  dz 


(<M  = 

1,2, 

6) 

(A-1.4) 

and  K 

ag  " 

:  Mg 

Mg 

;  (a,g  =  4,5) 

(A-1.5) 

where  [A1] 

= 

tL]"1 

H/2 

Mg 

= 

1 

H2 

/ 

Mgdz 

-H/2 

(A-1.6) 

[S] 

= 

[Q]-1 

tQ] 

_ 

I  Q44 

Q45 

m 

a 

°55_ 

and  k '  are 
ap 

the 

shear 

correction  factors  as 

described  in 

references 

24  ,  25. 
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APPENDIX  A- 2 

STRESS  ANALYSIS  IN  A  DISBONDED  LAMINATE 


In  this  section  a  method  of  solution  of  the  boundary  value 
problem  of  the  disbonded  laminate,  as  shown  in  figure  1,  is  out¬ 
lined.  Three-dimensional  elasticity  solution  is  not  attempted. 

A  lower  order  structural  theory  based  on  the  governing  differen¬ 
tial  equations  of  laminated  plates  outlined  in  Appendix  A-l  is 
utilized.  The  disbond  is  located  between  two  laminates  whose 
dimensions  in  x-  and  y-directions  are  infinitely  large.  The  sur¬ 
faces  of  the  disbond  may  be  subjected  to  surface  tractions  in  x-, 
y-,  and  z-directions .  Tractions  on  the  two  surfaces  are  equal  in 
magnitude  but  opposite  in  sign.  It  is  advantageous  at  this  stage 
to  introduce  the  following  non-dimensional  variables  correspond¬ 
ing  to  the  quantities  which  appear  in  the  differential  equations 
A-l. 3.  For  non-dimensionalization ,  use  is  made  of  two  quantities, 
i.e.  (i)  Hfl,  which  has  the  dimension  of  length,  and  (ii)  DQ,  which 
has  the  dimension  of  bending  rigidity. 


1.  Coordinates  -  [x1,x2l  =  [x,y]/Hp 


2  . 

Thickness 

-  h  =  H/Hq 

3. 

Displacement 

Variables  -  [U^U^U^] 

=  (u,v,w)/H0 

=  lW 

4. 

Stiffnesses 

-  *  aaSH0/4Do) 

Ke  ■  B«sVDoh  I0’8 

=  1,2,6 

daS  =  DaB^D0h  ' 

ea8  '  Vo'V  “'6 

=  4,5 

5. 

Prestress 

-  INU'f,22'N121  *  |NXX- 

N°  ,N°  ] H?/D. 
yy  xy  0  0 

6. 

Tractions 

-  [tx,  t2,  t3]  =  [t^ 

V  V  hJ/d0 

(A-2.1) 
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TRANSFORMED  VARIABLES 


Equations  (A-1.3)  can  now  be  written  in  terms  of  the  non- 
dimensionalized  variables  defined  in  (A-2.1)  and  reduced  to  five 
algebraic  equations  by  taking  multiple  Fourier  transforms. 
Transforms  of  variables  Uj ,  ^ ^ ,  and  t ^  are  denoted  with  a  bar  on 
top  of  the  quantity,  i.e. 


OO  CO 


Uj(sl'S2)  =  /  /  UjUl' 


i  (s1x1+s2x2) 


x2)  e 


dx^dx2  (A-2.2) 


-OO  — OO 


where  The  inverse  transform  is  given  by: 


Uj(xl'x2) 


OO  OO 

=  ^  /  /  vv 


-i  (s1x1+s2x2) 


s2)  e 


ds^ds2 . 


(A-2.3) 


Introduction  of  non-dimensionalized  displacement  variables  at 
z=+H/2  of  a  laminate  given  by: 


u  .  =  U  .  ±  .  h/2 ;  j=l ,  2 

3  1  3 


U3  =  U3 


(A-2.4) 


as  well  as  some  algebraic  manipulations  result  in  the  following 

_ + 

relationship  between  transformed  tractions  tT  at  z=±H/2  and  the 
displacements : 
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I 


1 

Hi? 

_ 1 

11 

12 

13 

1 - 

+  rH 

13 

1 

C 

C 

C 

fc2 

(2x2) 

(2x2) 

(2x1) 

U2 

-  — - 

—  — 

'q  1 

21 

22 

23 

"I 

C 

C 

C 

L  2  j 

(2x2) 

(2x2) 

(2x1) 

- 1 

1  CM 

’’i 

—  — 

.  _ 

31 

32 

33 

-i(t+  +  tZ) 

C 

C 

C 

-iu3 

(1x2) 

(1x2) 

(lxl) 

21 

Where,  [C  ] 


12  31 

[C  ]  and  [C  ] 


Tr 


32 

=  -1C  ] 


Tr 


23 
-1C  ) 


13 

[C  ]. 


ag 

The  elements  of  the  matrices  [C  ]  are  functions  of  laminate 
properties  as  well  as  s,  and  s  .  In  the  above  expression 
31  Tr 

[C  ]  indicates  transpose  of  [C  ] . 

For  the  purpose  of  carrying  out  the  algebraic  manipulations 
which  are  necessary  to  reduce  the  mixed  boundary  value  problem 
to  a  set  of  integral  equations,  the  following  transformations 
are  introduced. 


s^  =  t,cos<j>/£^ 
s2  =  csin<J)/f,2 

X1  =  {-^rcos0 
X2  =  ®-2rs^n® 


(A-2.6) 


(A-2.7) 
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where  and  are  two  parameters,  chosen  in  a  manner  such  that 
the  disbond  shape  can  be  mapped  on  to  a  region  which  is  close  to 
the  unit  circle  r<l.  In  particular  an  elliptic  shaped  disbond 
with  semiaxes  and  can  be  mapped  on  to  the  circle  r=l  by 
chosing 

aB 

In  terms  of  the  variables  c,4>  the  matrices  [C  ]  (a,  6=1, 2) 

can  be  written  in  the  following  fashion. 


aB 

fC  ] 


aB  aB 

2  1  2 

r[CXl  +  [<T] 


(A-2.8) 


where , 


22  12 

2  2 

[C  J  =  -[CT] 


55  S45 

45  e44 


(A-2.9) 


aB 

1 


'11 


aB  2  “B  2  »B 

Elnl  +  E2n2  +  2E3nin2 


'12 


aB  2  aB  2 

C21  =  E3ni  +  E4n2  +  E5nln2 


(A-2.10) 


aB 

1 

c 


22 


aB  2  aB  2  aB 

^2^1  +  E6r|2  +  ^^4r'lr,2 


n1  =  cos<J>/£^ 


(A-2.11) 


n2  =  sin<j>/£2 


-109- 


and 


11 

E1  =  all  +  ^dll^  +  ^bll^ 


11 

E2  =  a66  +  {d66}  + 


11 

E3  =  a16  +  ^d16  ^  +  (b165 


11 

E4  =  a26  +  {d26}  +  (b26) 


E5  ai2  +  a66  +  {d12+d66}  +  (b12+t>66) 


11 

E6  =  a22  +  d22  +  (b22) 


(A-2.12) 


22 

Expressions  for  E  ^  (j=l,6)  are  the  same  as  (A-2.12)  with  signs 
of  the  terms  within  parantheses  changed  to  negative.  To  obtain 
12  21 

Ej  =  Ej  (j=l,6)  it  is  necessary  to  drop  from  (A-2.12)  the  terms 
in  the  parantheses  and  change  the  signs  of  the  terms  in  brackets 
to  negative. 

33  13 

[C  ]  and  [C  ]  are  given  by: 


33 
[C  ] 


13 
(C  ] 


13 

ac1] 


(A-2.13) 
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33 

1  2 

where,  c  =  (e^+N^)^  +  2  (e45+N12^  )1in2  + 


13 

1 


'11 


(e55Ve45n2)/h 


13 

4 


(e45nl+e44n2^ //h 


TRACTION  FREE  CONDITION  ON  UPPER  SURFACE 


The  top  surface  of  laminate  1  (fig.  1) 
and  therefore: 


(e44+N22 


)  n 


2 

2 


(A-2.14) 


is  free  of  tractions 


(A- 2 .15) 


Elimination  of  u+  and  u+  from  (A-2.5)  results  in  an  expression 
of  the  following  form. 


[t-3 


— 

11 

12 

F 

F 

(2x2) 

(2x2) 

"2 

(i) 

L—  -J 

=  [  F  ] (u"l 

21 

22 

-  itT 

_ 

3 

F 

F 

-iu3 

(1x2) 

(lxl) 

i— .  — 

— 

_ 

_ 

(A-2.16) 
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where 


11 

22 

21  11  " 

1  12 

IF  ] 

=  [C  ]  - 

[C  ]  [C  ] 

[C  ] 

12 

21  Tr 

23 

21  11  _1  13 

[F  ] 

=  [F  ] 

=  [C  ]  - 

IC  ]  [C  ]  [C  ] 

22 

33 

31  11 

1  13 

[F  ] 

=  [C  ]  - 

[C  ]  [C  ] 

[C  ] 

The 

matrix  [F] 

can  also 

be  written  as 

(1) 

[  F 

(1) 

]  =  [  F*] 

♦  t'J*  , 

(A-2.17) 


(A-2.18) 


where  [F  *]  and  [F°]  are  of  the  following  form. 


(1) 

[  F*]  =  c' 


— ' 

— 

— 

— 

— 

fii 

4 

0 

0 

0 

f13 

fll 

f  12  ° 

f  21 

fl 

r22 

0 

+  Q 

0 

0 

f23 

+ 

f  21 

f  22  ° 

0 

0 

f33 

f  31 

f2 

32 

0 

0 

0  f  33 

— 

— 

— 

_ 

_  - 

(A-2.19) 


The  elements  f^D  (j  =  l,2,3)  are  functions  of  n,  and  r)->  and 

ap  12 
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[  F°] 


(A- 2 .20) 


f 

Cf 


21 

0 

31 


Cf 


22 

0 

32 


Cf 

f 


23 

0 

33 


The  elements 


f°.  can  be  written  as: 

a£3 


=  tt‘ 


+  <2‘*2  ♦  f*3’ 


(A- 2 .21) 


Oi  *  -j 

where  (j  =  l,2)  and  f  J  (j=l,2,3)  are  functions  of  and  n2 

only . 

LOWER  SURFACE  RESTING  ON  ELASTIC  FOUNDATION 

The  lower  surface  of  laminate  2  (fig.  1)  is  either  traction 
free  or  rests  on  an  elastic  foundation.  When  it  rests  on  elastic 
support  one  has  the  following  conditions. 


t~ 

3 


+  K°  uT 
3  3 


0  ;  j=l,2 


fc3  + 


K3  U3 


0 


(A- 2 .22) 


where  K ^  =  K^Hq/Dq,  (j=l,2,3)  being  modulus  of  the  elastic 
subgrade  (traction  to  cause  unit  deflection)  in  x- ,  y-  and 
z-directions ,  respectively.  Substitution  of  (A-2.22)  in  (A-2.5) 


and  elimination  of  u~  ar.d  u2  result  in  a  relationship  of  the  form 
given  below. 


(A- 2 .23) 


(2)  (1) 

where  [  F  ]  has  the  same  form  as  F  (see  equations  A-2.18  - 
A-  2.21) . 

CONDITIONS  AT  INTERFACE 

Transform  of  the  traction  vector  ( [t]=  [t^ ,  t2  ,-it-j] )  acting 

on  laminate  1  or  2  at  the  interface  containing  the  disbond  is 

related  to  the  transform  of  displacement  vector  ( [u]= [u^,u2,-iu3l ) 

of  the  same  laminate  at  the  same  interface  in  terms  of  the  stiff- 
(1)  (2) 

ness  matrices  [  F  ]  or  [  F  ]  by  means  of  equations  (A-2.16)  and 
(A-2.23).  In  what  follows  tractions  and  displacements  on 
laminates  1  and  2  will  be  identified  by  superscripts  -  and  +  , 
respectively.  Since  the  applied  tractions  on  top  and  bottom  sur¬ 
faces  of  the  disbond  (as  well  as  everywhere  else  on  the  interface) 
are  equal  in  magnitude  and  opposite  in  sign,  one  has  the  following 
identity . 


[t+]  = 

^2 

(2)  (2) 

=  [  F  ]  [u+]  =  [F] 

-it+ 

1 

"^3 

tf]  +  [t+]  =  0 


(A- 2 .24) 
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Furthermore,  if  (x^x^,  k=l,2,3  are  the  difference  between 
the  displacements  (at  the  interface)  of  laminates  1  and  2  in  x- , 
y- ,  z-directions ,  then: 


[u“]  -  [u+]  =  [u  ]  = 


[ul,u2,-iu3l 


(A-2.25) 


where  u£\,  k=l,2,3  are  the  transforms  of  u*.  The  four  equations 
(A-2.16)  and  (A-2 . 23-A-2 . 2 5)  can  now  be  formally  solved  and  the 
transformed  tractions  at  the  interface  can  be  expressed  as: 


[t“]  =  -[t+1  =  [G] [u*] 


(A-2. 26) 


where  [G]  is  given  by  the  following  relations. 


(1)  (1)  (2)  -1  m  „  -1  (2) 

[G  ]  =  [  F  ] [  F  +  F  )  [Fr]  +  [  F°] [F*]  [  F*] 


(1)  „  -1  (2) 

+  [  F*] [F* ]  [  F*] 


(i)  (i)  (i) 

[  F  ]  =  [  F*]  +  [  F°]  (as  in  (A-2. 18)) 


(A-2. 27) 


(1)  (2) 

[F*]  =  [  F*]  +  [  F*] 


(A-2. 28) 


and 


r  {2)n  ~  _1  (1) 

[Fr]  =  [  F°] IF*]  [  F*1 


(1)  _  -1  (2) 

[  F°] [F*]  (  F*] 


(A-2. 29) 
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The  matrix  [G]  can  also  be  expressed  as: 

[G]  =  C2  [G1]  +  UG21  +  [G3]  +  [G4]  (A- 2 .30) 

where  the  matrices  [G^]  (j=l,4)  are  symmetric  and  the  elements 

of  the  matrices  are  such  that: 


«i,  ■  4  -  4  ■  4  -  4  ■  ^12  ■  4  -  4  ■  »•  la-2-31’ 

Other  non-zero  elements  of  [G-1]  (j=l,2,3)  are  functions  of 

n1,n2  and  properties  of  the  two  laminates. 

4 

The  matrix  [G  ]  is  of  the  form: 


[G4]  = 


'll 


4 

*21 


5g 


31 


'12 


'22 


59 


32 


59 


59 


13 

4 

23 

4 

*33 


(A-2.32) 


where  g4g  =  g^g  (52) /g*  (52)  ,  <*,&  =  l,2,3;g4*  and  g*  being  poly- 

p  2 

nomials  of  5  with  coefficients  of  the  powers  of  c  ,  which  are 
functions  of  n1  and  n2  and  properties  of  the  laminates.  It  is 

*  4* 

noted  that  the  degree  of  the  polynomial  g  exceeds  that  of  gag 


by  unity . 
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A  SET  OF  INTEGRAL  EQUATIONS  FOR  ELLIPTIC  DISBOND 

The  displacements  of  the  two  laminate  :  the  interface 
differ  only  over  the  area  of  the  disbond  which  will  be  assumed 
to  be  an  ellipse  in  x^x2  plane  with  semiaxes  and  Further, 

the  ellipse  is  mapped  on  to  a  unit  circle  (r_<l)  by  the  trans¬ 
formation  (A-2.7).  It  is  possible  to  expand  the  discontinuity 
in  the  displacements  u£  across  the  disbond  as  well  as  the 
tractions  acting  on  the  top  surface  of  the  disbond  in  the  form: 


*  i 


w  2 

uk(r,0)  =  u*(x1,x2)  =  E  r  lpl  u*  (r  )eip6;  r  <  1;  k=l,2,3 

P  =  —  00  r 


=  0  r  >  1 


(A- 2 .33) 


00 

tT  (r,i|0  =  tT(x  ,x  )  =  E  t.  (r)eim*;  k=l,2,3  (A-2.34) 

3  3  -L  ^  m=-«> 


where  u£(_p)  and  tk.(-p)  are  comP^ex  conjugates  of  u*  ^  and 
t,  ,  respectively. 

k  p 

Taking  multiple  Fourier  transform  of  u^  (as  expressed  in 
equation  A-2.33)  results  in: 


00  oo 


u*(c,4>)  =  u*(s1,s2)  =  J  J u*  (x1,x2)e1  (slXl 


+  S2X2)dx1d^ 


—  00  —  oo 


1  2  TT 


-  I  l2  E  J  J  r 

n=-^  _ 


’ 2  \  QiPe„i^r  ' cos  (0-<j>) 


ukp(r  )e  ^  e 


0  0 


d6dr 

(A-2.35) 
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Use  of  the  identity  (ref.  26) 

2* 

r  eipeeur  cOS(6-«)d6,2lIe-ip*(i)PJ  (cr',  (A-2.36) 

**  P 

0 

yields 


nf  U  ,<t>) 


J, 

=  2-rrJ^Jl  I(i  )e~ip<t>f  rrlPl+1u*  (r2)J  (Cr’)dr' 

0=—  oo  .  b*  lr 


(A- 2 .37) 


where  J  denotes  Bessel  function  of  the  first  kind  of  order  p. 
P 

In  what  follows  the  following  relations  and  substitutions 
are  utilized. 


CO  CO 


(r  ,ijj)  =t~  (x  ,x  )  =  -Xr  f  [  tT(Sl  ,s,)e'l(slxl+s2x2)ds1ds, 
]  3  1  2  4it2  J  J  d  1  2  1  ^ 


—  CO  —  oo 


2tt 


-0  '0 


y  y  c  tj“  e  l!*rcos  d(^dc  (A-2.38) 


t.  (r) 
jm 


2n 

■s  / 1 


tj  (r  ,i(j)  e“im^di(j 


(A-2.39) 


2  TT 

/•  ^irtuj^-Urcos  <*-<(,)  dij,  =  27r(_A>™  e~^  j^r)  (a-2.40) 
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d  * 


,2 


,2 


)  =  r  u, _ (r  ) 


dr 


kp 


kp 


(A-2.41) 


Ukp(r’  >  =  Ukp  (r’2)  +  1  ukp  (r’2) 


(A-2.42) 


(r2>  +  1  <r2> 


(A-2.43) 


For  elements  g.^,  g12'  and  g33  matr;*-x  G  (in  equation 
A- 2 .30) 


2  IT  TT 

f  gj]c  ,  4> )  e1  (P“m)  ^d(J)  =  if  gjk  (c  ,♦>  e1  (p_m)  *d*  ;  p-m= 

**  r\  r\ 


even 


=  0 


p-m=odd 


(A- 2 .44) 


For  elements  g^  and  923  °f  matrix  [G] 


2tt 

f  gjk(C^)el(p“mUd<f  =  0 


;  p-m=even 


=  2  J  gjk  U  ,<j>)  e1  (p_m)  ^dcf) ;  p-m=odd 
0 

(A-2.45) 
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(A—  2.38)  and  (A-2.39)  are  the  laws  of  inversion  of  Fourier  trans¬ 
forms,  (A-2.40)  is  an  identity  (ref.  26)  and  (A-2.44)  as  well  as 
(A-2.45)  follow  from  the  properties  of  the  elements  of  of 
matrix  [G] .  Because  of  (A-2.44)  and  (A-2.45)  the  problem  may  be 
divided  into  two  parts  such  that  for  part  1: 


tR 

i  (2p) 

tl(1) 

lpl 

- 1 

rt 

M  M 

N) 

1 

t<2> 

lpl 

tI 

2(2p) 

= 

t(1) 

2P2 

9 

tR 

2(2p) 

= 

t(2) 

2P2 

tR 

t3(2p+l) 

a2p  Ul(2p) 

t(1) 

3p3 

u<1} 

lpl 

tI 

t3 ( 2p) 

!  ' 
a2p  Ul(2p) 

t(2) 

_  3p3_ 

u(2) 

lpl 

a2p  U2(2p) 

= 

u(1> 

2p2 

9 

2p  (2p) 

= 

u<2> 

2p2 

R 

U3(2p+1) 

3p3_ 

I 

U3(2p+1) 

U<2> 

-  3_ 

(A-2.46) 

and  for  part  2 : 


tR 

l(2p+l) 

I 

ft* 

H*  — 

TJ  t— 

J 

t1 

l(2p+l) 

t<2> 

lpl 

t1 

z2 ( 2p+l ) 

= 

t(1) 

t2p2 

9  “ 

tR 

2 (2p+l) 

t(2) 

2p, 

—  — 

-  — — 

—  — 

—  — 

R 

Ul(2p+1) 

“i(1) 

1P1 

i 

Ul(2p+1) 

u<2) 

lpl 

I 

U2(2p+1) 

= 

U^ 

2p2 

r 

R 

U2(2p+1) 

= 

(2) 
u'  ' 

2?  2 

a2p  U3(2p) 

u*1* 

_  P3  _ 

a2p  U3(2p) 

u<2> 

3p3 

(A-2.47) 


where  p  varies  from  0,1,2...°°,  and  a2p=l/2,  p=0;  a2p=l,  P>^- 
Substitution  of  (A-2.37)  in  (A-2.26)  as  well  as  use  of 
(A-2.38  -  A- 2 .47)  yield  the  following  set  of  equations  after 
integration  by  parts  with  respect  to  r'  and  utilization  of  the 
condition  u*p{l)=0. 


(ql} 
t.  (r) 


0  2  3 

-  Z  Z  Z  (-1) 

71  q2=l  k=l  p=0 


p,  -m . 
k  3 


(q,q 


1^*2 


/  Yjk(mj'V*)Jm.Ur) 


0  '0 


1 


/ 


r 


pkt2 


k+1 


(q2} 

u.  (r  2)dr  d$dc 
kpk 


(A-2.48) 


;  j  =  1,  . .  .3 
m  =  0,1,  ..  .0° 

qx  = 


where  the  repeated  indices  jk  do  not  indicate  summation  and 
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m . 
3 

=  2m  '  Pj  = 

2p 

for  j=l , 2 

;  for  part  1 

m . 

3 

=  2m+l ,  p .  = 

3 

2p+l  for  j=3 

m . 
3 

=  2m+lf  p.  = 

3 

2p+l  for  j=l,2 

;  for  part  2 

mj 

=  2m  '  Pj  = 

2p 

for  j  =  3 

(q^ 

q2) 

and  y 

ik  (ql'q2=1,2j 

j  . 

,k=l,2,3)  are  given  by  four  [3; 

as  given  below. 

(11) 

cosm^cosp^cj) 

-cosm^<J>sinp2iJ> 

cosm^cosp^ 

(Yl  = 

-sinm^cosp^ 

sinm2<{isinp2<J> 

-sinm2(j>cosp3(j) 

cosm3<j>cosp1  <f> 

-cosm3^sinp2<}> 

cosm3(f>cosp3<j> 

(12) 

-cosm3<Jisinp3<j> 

cosm3<J>cosp2<{> 

-cosm^sinp^ 

[y]  = 

sinm^sinp^tj) 

-sinm2<{icosp2(}> 

sinm2<)>sinp3<t> 

-cosm3<j)sinp^<t> 

cosm3$cosp24> 

-cosm3ij>sinp3<j> 

— 

(21) 

-sinm^cosp^ 

sinm^<))sinp2<!) 

-sinm^4>cosp3<J) 

(Yl  = 

cosm^cosp^ 

-cosm24)sinp2(J> 

cosm2({>cosp3<j> 

-sinm34>cosp^(}i 

sinm34>sinp24> 

-sinm34>cosp2(t> 

(22) 

sinm^sinp.^ 

-sinm^4)cosp2$ 

sinm^<J>sinp3<J> 

lY)  = 

-cosm24>sinp^(j) 

cosm^cosp^ 

-cosm2<))sinp3<J) 

sinm^sinp.^ 

-sinm^cosp^ 

sinm34)sinp3<{i 

(A- 2 .49) 


(A-2.50) 
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For  0  _<  r  <  1,  equations  (A-2.48)  are  a  set  of  integral 
equations,  which  must  be  solved  to  determine  the  displacements 
and  stress  states  in  the  disbonded  laminate. 


SOLUTION  OF  THE  SET  OF  INTEGRAL  EQUATIONS 


★ 

The  discontinuities  in  displacements  across  the  disbond, u^ 
can  be  calculated  with  the  help  of  equations  (A-2.46  or  A-2.47), 

(q2) 

(A-2.42),  (A-2.41),  and  (A-2.33)  after  the  functions  u,  are 

Kpk 


evaluated  by  solving  (A-2.48).  The  left-hand  side  of  (A-2.48) 
contains  functions  which  are  determined  from  the  prescribed 
tractions  on  the  disbond  surfaces  by  the  use  of  (A-2.39), 
(A-2.43),  and  (A-2.46  or  A-2.47).  To  solve  the  set  of  integral 


(q) 


o 

equations  (A-2.48),  the  functions  u,  (r1  )  are  expanded  in  a 

KPi, 


series  of  Legendre  polynomials  Pn_^(2r'  -1),  n=l,2,. 


(q)  2 


Z  u,  P  , 
n=1  kpqn  n-1 


(2r ' 2-l) 


(A- 2 .51) 


Substitution  of  (A-2 
tions  (A-2. 52  and  A-2. 54) 


51)  in  (A-2.48)  and  use  of  the  rela- 
given  below  (ref.  26): 


pn-l(2r'2-1> 


n  xn  n 


(1-r ' 2 ) 


"r1 


(A-2. 52) 
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Where 


'-In 


=  1  and  x-^  =  1 


'n^r. 


(-n+1)  ,(n)  i 

l"1  nl_1  * 

- n^TTi -  and  V 


=  (n1~l) ! ;  n1=l,2...n 


(n) 


n  _1  =  n  (n+1)  (n+2) . (n+n^- 


2) 


(A-2.53) 


Pu+2 


n,-l  n,-l  .  -n, 


r'  Jp  +1(?r') (1-r'2)  dr '=2  1  x„  ?  XJ 

k 


(O 


ni  pk+ni+1 

(A-2.54) 


as  well  as  integration  of  both  sides  of  the  resulting  expression 

m.+l 

between  0  to  r  after  multiplying  with  r  ^  results  in: 


(qx) 

t*  (r) 

3”j 


2  3 


OO  CO 


7  ^  2  £  z  (-1) 

q9=l  k=l  p=0  n=l 


Pk“mj 


Vn^lS"2  "  r 


n.,-1  nj.+l 
1  -  3 


CO  J 

// 

0  o 


-ni-i 

93k(C,<f>)  Yjk(n’j,pk'<t’)  Jpk+ni  +  l(C)  Jm;.  +  l(;r)d<))di; 


(A-2.55) 
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i 


Where : 


(q)  r  (q)  m.+l 

t*  (r)  =  /  t  .  (r')r'  3  dr1 

]m.  J  Dm. 

J  0  3 


A-2.56) 


If  one  expands  t*m_.  (r)  ;  0  <  r  <  1  in  a  series  of  Legendre  poly- 
2  ' 

nomials  Pn,_^(2r  -1)  ,  n  =1,2.  .  .<=°,  then  the  coefficients  are 
given : 


1  (q, ) 


1  a 

SmV  '  4J  r  ‘‘jm.  <r) 


2  3  °°  00  j  mq  1  n ' 

Z  Z  Z  Z  ,  1  u. 

kpq?n  kpq-,n 

,=1  k=l  p=0  n=l  2  2 


<  j  -1 »  2 , 3 
m  =0,1,2. 
q1=l,2 


n  '  =  1 , 2  , 


(A-2.57) 


Equation  (A-2.57)  is  an  infinite  system  of  linear  algebraic 

jmq,n’ 

equations  involving  the  unknowns  u,  .  The  coefficients  , , 

kpq2n  Lkpq2n 

are  given  by: 


jmq  n'  p.-m  n  n' 

x  =  -  (-1)  K  3  T.  r. 


jkpq2n  tt 


1  1  Xn  nxn  ’n '  2 

n1=l  n^=l  nln  nln 


nl+ni 


(A-2.58) 


/  /  1 


n.-n'-i 


(qiq2} 


0  0 


gjkU,l),)  Yjk  {mj'Pk'<!,)Jpk+n1+lU)Jmj+n^+l(c)d<})dc 


g.,  are  the  elements  of  the  matrix  [G]  which  is  given  in 
1  ^ 

equation  (A-2.31)  and  can  be  rewritten  as  the  sum  of  four  terms, 
i.e.  : 


[G]  =  Z  C3"6  [Ge((j>)]  +  [G4U,4>H 
0=1 


(A-2.59) 


Substitution  of  (A-2.59)  in  (A-2.58)  results  in  four  infinite 
integrals,  the  first  three  of  which  can  be  evaluated  in  closed 
form  by  the  use  of  the  following  relation  (ref.  26  ) : 


r  ,  rmru.) 

J  xT  JjOJjOd?  =  i 


U  V 


2Ar(x2)r(x3)r(x4) 


7  X^ ,  X^  ^  ^ 


=0  ;  X2  or  X4  £  0  (A-2.60) 

where 

X^  =  (v+u-X+l)/2 
X2  =  (-v+y+X+1) /2 
X^  =  ( u+ p+ X+l ) / 2 
X^  =  (v-p+X+D/2 

and  T  denotes  the  Gamma  function.  It  is  necessary  to  evaluate 
the  integral  (A-2.60) ,  only  when  X  and  X^;  8=1,... 4  are  positive 
integers.  Integration  with  respect  to  $  can  be  performed  numer¬ 
ically  by  the  formula  given  below. 


0(*)d*  =  rp  £  ut  fi(<j>R  ) 
Mi  e^o  6i  6i 


(A- 2 .61) 


where 


Ug  «  1/2  ,  0L  =  0fMx 


,  $  =  1,2. . .M1-l 


♦  6  =  ir01/M1  ,  6X  =  0,1 . Ml. 

4 

The  last  integral,  which  is  the  contribution  due  to  [G  ]  in 
(A-2.59  and  A-2.30)  contains  more  complicated  functions  of  c 

4 

and  4>  •  In  this  study  [G  ]  is  evaluated  numerically  by  the  use  of 
(A-2.27)  for  various  values  of  5  and  <f>  and  the  integral  is  eval¬ 
uated  numerically  by  the  formula: 


71 

.  /  / 


AU  ,<)>)d<j>d(;  = 


£  E  I  w?  A(Cr  ) 
M1B2=1  B1=0  32  B1  e2  B1 


(A-2.62) 


where , 


1-5, 


2 

Cg  =  - 5( 

P2  1  +  5d 


•s2  -K/,i+v2:k 


CR  are  the  zeroes  of  Legendre  polynomic  PM  (C  ) ;  -1  <  C  <1 
^2  2 

* 

and  w  are  the  weights  for  Gaussion  integration  and  is  a  multi- 
°2  0 
plier  chosen  for  obtaining  rapid  convergence  of  the  integrals. 
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For  the  purpose  of  numerical  solution  it  is  necessary  to 
truncate  the  infinite  system  (A-2.57): 


m. 


p  <  m 
*  —  max 


n'  n’  1  "max' 


Further,  for  laminates  which  are  orthotropic  with  x,y  axes  as 
the  principal  material  axes,  the  set  of  equations  for  ukpln  and 
ukp2n  are  uncouPled  and  can  be  solved  separately. 

CALCULATIONS  OF  STRAIN  ENERGY  RELEASE  RATES 

It  follows  from  equations  (A-2.33),  (A-2.41),  (A-2.42), 
(A-2.46),  (A-2.47)  and  (A-2.51)  that  the  derivatives  of  u*(r;6) 
as  r'  approaches  1  are  given  by: 


3u*  |  00  00 

IfH  ,  =  2  E  En  [ukplnCOSPk9'Ukp2nSinpk9] ;  k=1'3 

r  +1  p=0  n=l 


oo  oo 


2  E  [-ukplnsinPke+Ukp2ncosPk01;  k=2 

p=0  n=l  ^ 


(A- 2. 


Denoting  the  displacement  discontinuities  in  directions  normal 
and  tangential  to  the  elliptical  boundary  of  the  disbond  by  Uj 
and  Ujj,  respectively,  and  in  z-direction  by  UII:I»  the  rate  of 
change  of  these  quantities  at  the  boundary  in  the  direction  x^, 
normal  to  the  boundary  curve  are  given  by: 


63) 
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(A- 2 .64) 


r, ,  _  SIIIO  j.  tosd  z 

II  “  sl2  3r’  3r' 


Uill 


=  / 


'  2 
cos  0 


- 2~ 

sin  0 

~T~ 


3u* 

‘3rT 


Since  these  derivatives  are  zero  outside  the  disbond  area,  they 
have  a  discontinuity  as  the  boundary  is  crossed.  As  will  be  seen 
next,  the  integrals  of  tractions  imposed  by  one  of  the  laminates 
on  the  other  at  the  plane  of  the  disbond  with  respect  to  the  nor¬ 
mal  Xj  also  has  a  discontinuity,  indicating  existence  of  concen¬ 
trated  (line)  forces  at  the  disbond  boundary.  This  is  the  type 
of  force  singularity  which  the  laminated  plate  theory  assumptions 
yield  in  contrast  to  the  type  of  stress  singularity  encountered 
in  2-D  or  3-D  elasticity  theories. 

Substitution  of  (A-2 . 51-A-2 . 54)  and  (A-2.30)  in  (A-2.48)  and 
omission  of  all  terms  which  do  not  yield  concentrated  line  forces 
on  the  elliptic  boundary  (r=l)  yields  the  following  expression 
from  which  the  magnitudes  of  the  concentrated  forces  are  to  be 
obtained . 


(r) 


L  L 

k  p=0 


E 

n=l 


pk"mi 

(-1)  :u 


kpq2n 


00 


0 


/  (^lc^2) 

j  gjk(*}  Yjk(mj'pk'^ 


Q  J 


.,U)  J„  (Cr)  d$dC 


(A-2. 65) 


For  j=l,2  the  summation  over  k  is  to  be  performed  over  k=l,2  and 
for  j=3  k  is  taken  equal  to  3.  This  follows  from  (A-2.31). 
Integration  with  respect  to  $  can  be  performed  numerically  by  the 
formula  (A-2.61).  The  infinite  integral  (w.r.t.  £)  can  be  written 


/  «  V+2('-IJm.<Sr,dt 
0  k  3 


m .  -1 
3  d 

-  r  [ft(pk,itK  »r)  ] 


(A-2.66) 


where 


ft(pk,m..,r)  =  r  J  Jp  +2U)  J^SrJdi;  ;  itk  =  0  (A-2.67) 


H .  r 

3  I  J  ,9U)  j  .(cr)dc;  m  >0  (A-2.68) 

J  p,+2  m.-l  3 


The  expressions  (A-2.67)  and  (A-2.68)  are  discontinuous  at 
r=l  and  can  be  expressed  as  (ref.  26) : 


I 


fi(pk,0,lT)  =  0 


fi(pk,0,l  ) 


Pk+4 

HzrL  F* 

p  +2  (  2  ' 

r  (-~) 


~  ,  2,  1) 


and  for  m.>0 
D 


p  +m . +2 

- )  Pk+m.+2  m.-p  -4 

a(pk,m.,l+)  =  -  -„-p  -  pbl-S-J— .  - 

T  i  3  \  ■p  /  — 


ir—‘ '  pk+3'^ 


r  ( — 1 — — )  r  (pk+3) 


;  nij  >  Pk  +  2 


=  0  ;  n»j  1  pk  +  2 


Pk+m-i+2 
r(^-J — ) 


fi(pk/  Wj  ,1  )  -  p  in, +4 


,  p,  +m  .  +  2 

Fh(^V- 


pk-m  +2 


-)r(m.) 


;  m .  <  p,  +4 


=  0  ;  m .  >  p,  +4 
3  -  *k 


(A- 2 .69) 


where  Fn  ,u2 ,y2 ,z)  is  the  hypergeometric  series  with  y^,  y2 
and  y^  as  integers,  which  terminates  because  of  the  fact  that  y2 
is  negative. 

If  Xj  is  the  direction  normal  to  the  elliptic  boundary  with  i 

semiaxes  and  then  the  partial  derivative  of  a  function  A 
with  respect  to  x^.  is  given  by:  / 


2 

sin  6 


2  2 

sin6cos0  ( i^+SL2)  ^  ^ 

n — 2 - 2 — r 

l  1 1 2c°s  0  +  Ji^sin  0 


(A-2.70) 


Therefore,  the  concentrated  force  components  in  nondimens ional 
form  are  given  by : 


(qx) 

t.  ' 

3m 


2  QO  OO 

E  Z  E  E  (-1) 


Pk-mj 


j  /  2 . 2Z  q  =1  k  p=0  n=l 

J  „/cos  0  .  sin  0  ^2  r 

TT/  - s—  +  - o— 


kpq2-n 


[fi(pk,mj  ,l+)-n(pk,m:j  ,1")  ]  /  g  (<J>)  yjk  (pk,mj  ,  4» )  d<(> 

0 


i 

(A- 2 .71) 


i 


1 

1 
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Summing  over  m  and  q^,  one  has  the  following  expression  for 
forces  in  x-,  y-,  and  z-directions : 


t  =  £  2a  [-t.  cosm.6  +  t.  sinm.0]  ;  j  =  l,3 


m=0  j 


m.  3m 


=  Z  2am .  (tjm.  sinmj6  “  fcjm.  cosmj03  ;  j-2 
m=0  3  3  3 


where,  a  =  1/2;  m.  =  0 

mj  3 


=  1  ;  m..  >  0 


(A- 2 .72) 


If  Tj  and  Tjj  are  the  forces  in  directions  normal  and  tan¬ 
gential  to  the  elliptic  boundary,  respectively,  and  is  the 
force  in  z-direction  (in  units  of  force) ,  then: 


,cos0  tc  ,  sin0 
'  ^  C1  z2' 


sin0  ^ c  ,  cos0  ,cx 


(A- 2 .73) 


/  27  T27 

,  .  _  /  cos  0  ,  sin  0 

where  A  =  ✓ - ^ 


and  Dq  and  are  the  quantities  introduced  for  nondimensionali- 
zation  (see  equation  A-2.1). 


The  strain  energy  release  rates  for  modes  I,  II,  and  III  can 
then  be  calculated  at  various  positions  on  the  elliptic  boundary 
by  the  relationship: 


